FROM BRAKE TO SYZYGY 
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Abstract. In the planar three-body problem, we study solutions with zero 
initial velocity (brake orbits). Following such a solution until the three masses 
become coUinear (syzygy), we obtain a continuous, flow-induced Poincare map. 
We study the image of the map in the set of coUinear configurations and define 
a continuous extension to the Lagrange triple collision orbit. In addition we 
provide a variational characterization of some of the resulting brake-to-syzygy 
orbits and find simple examples of periodic brake orbits. 



1. Introduction and Main Results. 

This paper concerns the interplay between brake orbits and syzygies in the New- 
tonian three body problem. A brake orbit is a solution, not necessarily periodic, for 
which the velocities of all three bodies are zero at some instant, the 'brake instant'. 
Brake orbits have zero angular momentum and negative energy. A syzygy occurs 
when the three bodies become coUinear. We will count binary collisions as syzygies, 
but exclude triple collision. 

Lagrange [5] discovered a brake orbit which ends in triple collision. The three 
bodies form an equilateral triangle at each instant, the triangle shrinking homo- 
thetically to triple collision. Extended over its maximum interval of existence, 
Lagrange's solution explodes out of triple collision, reaches maximum size at the 
brake instant, and then shrinks back to triple collision. 

Lagrange's solution is the only negative energy, zero angular momentum solution 
without syzygies |16lll7j . In particular, brake orbits have negative energy, and zero 
angular momentum, and so all of them, except Lagrange's, suffer syzygies. Thus 
we have a map taking a brake initial condition to its first syzygy. We call this 
map the syzygy map. Upon fixing the energy and reducing by symmetries, the 
domain and range of the map are topologically punctured open discs, the punctures 
corresponding to the Lagrange orbit. Sec figure [4] where the map takes the top 
"zero- velocity surface", or upper Hill boundary, of the solid Hill's region to the 
plane inside representing the coUinear configurations. The exceptional Lagrange 
orbit runs from the central point on the top surface to the origin in the plane 
(which corresponds to triple collision), connecting the puncture in the domain to 
the puncture in the range. 

Theorem 1. The syzygy map is continuous. Its image contains a neighborhood 
of the binary collision locus. For a large open set of mass parameters, including 
equal masses, the map extends continuously to the puncture, taking the equilateral 
triangle of Lagrange to triple collision. 
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Remark on the Range. Numerical evidence suggests that the syzygy map 
is not onto (see figure [5]). The closure of its range lies strictly inside the coUinear 
Hill's region. A heuristic explanation for this is as follows. The boundary of the 
domain of the syzygy map is the coUinear zero velocity curve, i.e., the coUinear Hill 
boundary. Orbits starting on this curve remain coUinear for all time and so are 
in a permanent state of syzygy. Nearby, non-coUinear orbits oscillate around the 
coUinear invariant manifold and take a certain time to reach syzygy, which need not 
approach zero as the initial point approaches the boundary of the domain. Thus the 
nearby orbits have time to move away from the boundary before reaching syzygy. 
It may be that the syzygy map extends continuously to the boundary but we do 
not pursue this question here. 

Collision- free syzygies come in three types, 1, 2, and 3 depending on which 
mass lies between the other two at syzygy. Listing the syzygy types in temporal 
order yields the syzygy sequence of a solution. (The syzygy sequence of a periodic 
collision-free solution encodes its free homotopy type, or braid type.) In ([H], 
the notion of syzygy sequence was used as a topological sorting tool for the three- 
body problem. (See also [19] and [2].) A "stuttering orbit" is a solution whose 
syzygy sequence has a stutter, meaning that the same symbol occurs twice in a 
row, as in "11" "22" or '33". For topological and variational reasons, one of us 
had believed that stuttering sequences were rare. The theorem easily proves the 
contrary to be true. 

Corollary 1. Within the negative energy, zero angular momentum phase space 
for the three body problem there is an open and unbounded set corresponding to 
stuttering orbits. 

Proof of the corollary. If a coUinear configuration q is in the image of the syzygy 
map, and if v is the velocity of the brake orbit segment at g, then by running this 
orbit backwards, which is to say, considering the solution with initial condition 
(q, — w), we obtain a brake orbit whose next syzygy is with velocity -\-v. This 
brake orbit is a stuttering orbit as long as q is not a collision point. Perturbing 
initial conditions slightly cannot destroy stutters, due to transversality of the orbit 
with the syzygy plane. QED 

Periodic Brake Orbits. In 1893 a mathematician named Meissel conjectured 
that if masses in the ratio 3, 4, 5 are place at the vertices of a 3-4-5 triangle and 
let go from rest then the corresponding brake orbit is periodic. Burrau [3] reported 
the conversation with Meissel and performed a pen-and-paper numerical tour-de- 
force which suggested the conjecture may be false. This "Pythagorean three-body 
problem" became a test case for numerical integration methods. Szehebely [27] 
carried the integration further and found the motion ends (and begins) in an elliptic- 
hyperbolic escape Peters and Szehebeley [28] perturbed away from the Pythagorean 
initial conditions and with the help of Newton iteration found a periodic brake orbit. 

Modern investigations into periodic brake orbits in general Hamiltonian systems 
began with Seifert's |22j 1948 topological existence proof for the existence of such 
orbits for harmonic-oscillator type potentials. (Otto Raul Ruiz coined the term 
"brake orbits" in [21].) We will establish existence of periodic brake orbits in 
the three-body problem by looking for brake orbits which hit the syzygy plane C 
orthogonally. Reflecting such an orbit yields a periodic brake orbit. Assume the 
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masses are nii = m2 = l,m3 > 0. Then there is an invariant isosceles subsystem 
of the three-body problem and we will prove: 

Theorem 2. For in an open set of mass parameters, including — 1, there 
is a periodic isosceles brake orbit which hits the syzygy plane C orthogonally upon 
its 2nd hit (see figure^. 

Do their exist brake orbits, besides Lagrange's whose 1st intersection with C is 
orthogonal? We conjecture not. Let I{t) be the total moment of inertia of the three 
bodies at time t. The metric on shape space is such that away from triple collision, 
a curve orthogonal to C must have J = at intersection. This non-existence 
conjecture would then follow from the validity of 

Conjecture 1. I{t) < holds along any brake orbit segment, from the brake time 
up to and including the time of 1st syzygy. 

Our evidence for conjecture [l] is primarily numerical. If this conjecture is true 
then we can eliminate the restriction on the masses in theorem [ij See the remark 
following proposition [Tol 

Variational Methods. We are interested in the interplay between variational 
methods, brake orbits, and syzygies. If the energy is fixed to be —h then the natural 
(and oldest) variational principle to use is that often called the Jacobi-Maupertuis 
action principle, described below in section |4] (See also [2], p. 37 eq. (2).) The 
associated action functional will be denoted Ajm (see eq. ([BO])). Non-collision 
critical points 7 for Ajm which lie in the Hill region are solutions to Newton's 
equations with energy —h. Curves inside the Hill region which minimize Ajm 
among all compenting curves in the Hill region connecting two fixed points, or two 
fixed subsets will be called JM minimizers. 

We gain understanding of the syzygy map by considering JM-minimizers con- 
necting a fixed syzygy configuration q to the Hill boundary. 

Theorem 3. (i).JM minimizers exist from any chosen point qo in the interior of 
the Hill boundary to the Hill boundary. These minimizers are solutions. When not 
collinear, a minimizer has at most one syzygy: go- 

(ii) . There exists a neighborhood lA of the binary collision locus, such that if qa € 
U, then the minimizers are not collinear. 

(iii) .IfqQ is triple collision then the minimizer is unique up to reflection and is 
one half of the Lagrange homothetic brake solution. 

Proof of the part of Theorem [l] regarding the image. Let W be the 

neighborhood of collision locus given by Theorem 3. If qo (z U, the minimizers of 
Theorem 3 realize non-collinear brake orbits whose first syzygy is qQ, therefore the 
image of the syzygy map contains U. QED 

An important step in the proof of theorem [3] is of independent interest. 

Lemma 1. [Jacobi-Maupertuis Marchal's lemma] Given two points qo and qi in 
the Hill region, a JM minimizer exists for the fixed endpoint problem of minimizing 
Ajm{i) among all paths 7 lying in the Hill region and connecting qo to qi. Any such 
minimizer is collision-free except possibly at its endpoint. If a minimizer does not 
touch the Hill boundary (except possibly at one endpoint) then after reparametriza- 
tion it is a solution with energy —h. 
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We can be more precise about minimizers to binary collision when two or all 
masses are equal. Let rij denote the distance beween mass i and mass j. 

Theorem 4 (Case of equal masses.), (a) Ifnii = m2 and if the starting point qq is 
a collision point with ri2 = then the minimizers of Theorem^ are isosceles brake 
orbits: ri3 = r23 throughout the orbit. 

(b) If mi — m2 and if the starting collinear point qq is such that ria < r23 (resp. 
I'lS > T2'i), then a minimizerj of Theorem^satisfies this same inequality: at every 
point "f{t) we have ri3(t) < r23{t) (resp. ri^it) > r23,{t)). 

(c) If all three masses are equal, and if qo is a collinear point, a minimizer 7 of 
Theorem satisfies the same side length inequalities as qo : if ri2 < ^13 < r23 for 
qo, then at every point j{t) of j we have ri2(t) < ri3{t) < r23(t). 

Part (c) of this theorem suggest: 

Conjecture 2. // three equal masses are let go at rest, in the shape of a scalene 
triangle with side lengths ri2 < ris < r23 and attract each other according to 
Newton's law then these side length inequalities f'i2(i) < ^13 (0 < ^23 (i) persist up 
to the instant t of first syzygy, 

Commentary. Our original goal in using variational methods was to construct 
the inverse of the syzygy map using JM minimizers. This approach was thwarted 
due to our inability to exclude or deal with caustics: brake orbits which cross each 
other in configuration space before syzygy. Points on the boundary of the image of 
the syzygy map appear to be conjugate points - points where non-collinear brake 
orbits "focus" onto a point of a collinear brake orbit. 

Outline and notation. In the next section we derive the equations of motion 



in terms suitable for our purposes. In section 3.1 we use these equations to rederive 
the theorem of [16], [17] regarding infinitely many syzygies. We also set up the 
syzygy map. In section [3. 3| we prove theorem 1 regarding continuity of the syzygy 
map. In section [4] we investigate variational properties of the Jacobi-Maupertuis 
metric and prove theorems [3] and |4] and the lemmas around them. In section [5] we 
establish Theorem 2 concerning a periodic isosceles brake orbit. 

2. Equation of Motion and Reduction 

Consider the planar three-body problem with masses mi > 0, i = 1, 2, 3. Let the 
positions be qi £ IR? = C and the velocities be Vi — qi G M^. Newton's laws of 
motion are the Euler-Lagrange equation of the Lagrangian 

(1) L^K + U 
where 

2K = mil^ip + TO2IW2P + ^-3^3^ 

(2) ^ _ mim,2 ^ mim^ ^ ma 7713 

ri2 ri3 r23 

Here = \qi — qj\ denotes the distance between the i-th and j-th masses. The 
total energy of the system is constant: 

K-U^-h h>0. 

Assume without loss of generality that total momentum is zero and that the 
center of mass is at the origin, i.e., 

miVi + m2V2 + m3V3 = miqi + 771252 + "^293 = 0. 



BRAKE-SYZYGY 



5 



Introduce Jacobi variables 

niiqi + m2(72 



(3) Ci = '?2 - qi 6 = 93 



mi + m2 



and their velocities ^i. Then the equations of motion are given by a Lagrangian of 
the same form ([I]) where now 

(4) ^ _ mim2 ^ mims ^ 17121713 

ri2 ri3 r23 

The mass parameters are: 

mim2 (mi+m2)m3 (mi + m2)m3 

(5) fii = /i2 = — 

mi + m2 mi + m2 + ms m 

where 

m = mi + m2 + 
is the total mass. The mutual distances are given by 

ri2 = 161 

(6) ri3 = IC2 + 1^26! 

'"23 = 16 - I'lCll 

where 

mi m2 

Pl = V2 



Till +m2 mi + m2 

2.1. Reduction. Jacobi coordinates ^ eliminate the translational symmetry, re- 
ducing the number of degrees of freedom from 6 to 4. The next step is the elimi- 
nation of the rotational symmetry to reduce from 4 to 3 degrees of freedom. This 
reduction is accomplished by fixing the angular momentum and working in the quo- 
tient space by rotations. When the angular momentum is zero, there is a particulary 
elegant way to accomplish this reduction. 

Regard the Jacobi variables ^ as complex numbers: ^ = (6j6) ^ C^. Introduce 
a Hermitian metric on C^: 

(7) ((u, v)) = fiiviwi + H2V2W2 

If = {{v,v)) denotes the corresponding norm then the kinetic energy is given 

by 

ll^lP = ml6P + M2l6P = ^ 

while 

(8) |ieip = Mii6p + /^2i6P = / 

is the moment of inertia. We will also use the alternative formula of Lagrange: 

(9) ll^lp = ^ (mim2ri2 + TniTTisrj^ + m2m3r23) 

where the distances r.y- are given by ([6]). The real part of this Hermitian metric 
is a Riemannian metric on C^. The imaginary part of the Hermitian metric is 
a nondegenerate two-form on with respect to which the angular momentum 
constant to of the three-body problem takes the form 

^ = im((60)- 
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The rotation group = 50(2) acts on according to (^1,^2) 6*^(6,^2)- 
To ehminate this symmetry introduce a new variable 

to measure the overall size of the configuration and let [£] = [^1,^2] G CP^ be 
the point in projective space with homogeneous coordinates ^. Explicity, [f] is an 
equivalence class of pairs of point of \ where ^ = ^' if and only if ^' — for 
some nonzero complex constant k. Thus [£] describes the shape of the configuration 
up to rotation and scaling. The variables (r, [£]) together coordinatize the quotient 
space CV^^ 

Recall that the one-dimensional complex projective space CP^ is essentially 
the usual Riemann sphere C U cx). The formula a ([^1,^2]) = C2/C1 gives a map 
a : CP^ — > C U CXI, the standard "afhne chart". Alternatively, one has the diffeo- 



morphism St : CP^ — >■ S"^ (defined in (33 1 ) to the standard unit sphere S"^ C Ffi by 
composing a with the inverse of a stereographic projection map cr : S*^ — )• C U 00. 
The space S = CP^ in any of these three forms will be called the shape sphere. 

In the papers [16] and [5] the sphere version of shape space was used, and the 
variables r, [^] were combined at times to give an isomorphism C'^/S'^ send- 
ing (r, [^]) 1-^ rSt{[^]). The projective version of the shape sphere, although less 
familiar, makes some of the computations below much simpler. Triple collision cor- 
responds to e and the quotient map \ — ?• (C^ \ Q)/S^ is realized by the 
map 

(10) ^ : C2 \ ^ Q (0,^) X Cpi;7r(0 = {Ml K]). 

To write down the quotient dynamics we need a description of the kinetic energy 
in quotient variables, and so we need a way of describing tangent vectors to CP^. 
Define the equivalence relation = by (^, u) = (^', it') if and only if there are complex 
numbers fc, I with fc 7^ such that (^', u') = (A:^, ku + It is easy to see that two 
pairs are equivalent if and only if T7r(^,w) = Ttt(£^\u') where Tn : T(C^ \ 0) ^■ 
TCP^ is the derivative of the quotient map. Thus an equivalence class u] of such 
pairs represents an element of T^^jCP , i.e., a shape velocity at the shape [^]. One 
verifies that the expression 

(11) iiK,c]iP = ||Ulia6-6ar 

defines a quadratic form on tangent vectors at [^] and as such is a metric. It is 
the Fubini-Study metric, which corresponds under the diffeomorphism St to the 
standard 'round' metric on the sphere, scaled so that the radius of the sphere is 
1/2. We emphasize that in this expression, and in the subsequent ones involving 
the variables r, [^], the variable f is to be viewed as a homogeneous coordinate on 
CPI so that the lieil = (Miieil'+M2|6l')'^' 

occurring in the denominator is not 



linked to r, which is taken as an independent variable. Indeed (11) is invariant 
under rotation and scaling of ^,1^. 

We have the following nice formula for the kinetic energy: 



1 ( „-,2 , ^ , „2 1 



Proposition 1. = = | | r" H — ^ + r" 
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We leave the proof up to the reader, or refer to [5] for an equivalent version. 
Taking uj = gives the simple formula 

(12) Ko^lf' + y'^\U2~aii\' 

By homogeneity, the negative potential energy U can also be expressed in terms 
of r, [C]. Set 

(13) n[e]) = iieiit^(0- 

thus defining V. Equivalently, X^([C]) = t^(C/ll'?ll)- Since the right-hand side is 
homogeneous of degree with respect to ^, and since U is invariant under rotations, 
the value of T^([C]) is independent of the choice of representative for [£]. Clearly we 
have U{0 = ^^(Kl) for f e \ and (r, [^]) = 7r(C). The function V : S ^ M 
will be called the shape potential. 

The function Lred ■ TQ — > IR given by 

(14) L,ed(r,f,K,e])=i^o + -V'(K]) 

r 

will be called the reduced Lagrangian. The theory of Lagrangian reduction [TTl [T5] 
then gives 

Proposition 2. Let ^(t) be a zero angular momentum solution of the three-body 
problem in Jacobi coordinates. Then {r(t), [^{t)]) = 7r(^(t)) € Q = (0,oo) x CP is 
a solution of the Euler- Lagrange equations for the reduced Lagrangian L^ed on TQ. 

2.2. The Shape Sphere and the Shape Potential. Let C be the set of coUinear 
configurations. In terms of Jacobi variables = (Ci> C2) G C if and only if the ratio 
of £,1,5,2 is real. The corresponding projective point then satisfies [£,] € RP^ C CP"'^. 
Recall that RP^ can also be viewed as the extended real line MUoo or as the circle 

. Thus one can say that the normalized coUinear shapes form a circle in the 
shape sphere, S. Taking the size into account one has C = x C x S. 
Here and throughout, by abuse of notation, we will write C as the set of coUinear 
states, either viewed before or after reduction by the circle action (so that C C Q), 
or by reduction by the circle action and scaling (so C C CP^). 

The binary collision configurations and the Lagrangian equilateral triangles play 
an important role in this paper. Viewed in CP^, these form five distinguished 
points, points whose homogeneous coordinates are easily found. Setting the mutual 
distances ^ equal to zero one finds collision shapes: 

&12 = [0,1] &13 = [1,-^^2] h23^[l,Vi] 

where, as usual, the notation [^1,^2] means that (^1,^2) is a representative of the 
projective point. Switching to the Riemann sphere model by setting z = ^2/^1 gives 

612 = 00 fel3 = -V2 b23 = 1^1- 
The equilateral triangles are found to be at [1, G CP^ or at l± £ C where 

, mi-m2 , \/3. V1-V2 , Vs. 

(15) l± = r ± 1^ * = ^ ± 1^ 

^ ' 2(mi+m2) 2 2 2 

We will choose coordinates on the shape sphere such that all of these special 
shapes have simple coordinate representations 16^ In these coordinates, the shape 
potential will also have a relatively simple form. To carry out this coordinate 
change, we use the well-known fact from complex analysis that there is a unique 
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conformal isomorphism (i.e., a fractional linear map) of the Riemann sphere taking 
any triple of points to any other triple. Thus one can move the binary collisions 
to any convenient locations. We move them to the third roots unity on the unit 
circle. Working projectively in homogeneous coordinates, a fractional linear map 

cw + d 



becomes a linear map 



where = [1,^] and [771,772] = 

Proposition 3. Let X — and 



"6' 




a b 










c d 




m 



he the unique conformal map taking 1, A, A to 
^12,^13,^23 respectively. Then cj) maps the unit circle to the collinear shapes and 
and 00 to the equilateral shapes 1+ and Z_ . Moreover, in homogeneous coordinates 



(16) 



= 



" 1 -l' 




m 


1+ 







ly. 



The proof is routine, aided by the fact that preserves cross ratios. 

We will be using 77 = (771,772) as homogeneous coordinates on CP^ and setting 

w = 772/771 e C. 

In w coordinates, we have seen that the collinear shapes form the unit circle, with 
the binary collisions at the third roots of unity and the Lagrange shapes &i w — 
0, 00. We also need the shape potential in variables. It can be calculated from 
the formulas in the previous subsection by simply setting 

(17) = m - m 6 = '+7/1 - L772 ?7i = 1 

First, ^ gives the remarkably simple expressions 

•2z-{x+\f + {y- 

Using these, one can express the norm of the homogeneous coordinates and the 
shape potential as functions of {x,y). ||^|| is given by ^ and 

f 77117772 , 771,1777,3 7772 W3 A 7711777,2 7)1,17773 7)727^13 



(18) 



' 13 
„2 



y3N2 
2 > ■ 



(19) V{x,y) = m\ 
with p,j = 



V ^12 



^■13 



^"23 / 



P12 



Pis 



P23 



Remark. It is worth saying a bit about the meaning of the expressions eq. ( 18 ) 



and the variables pij occuring in eq (19). A function on which is homogeneous 
of degree and rotationally invariant defines a function on CP^. But the are 
homogeneous of degree 1, so do not define functions on CP^ in this simple manner 
So what is eq. ( 18 ) saying? Introduce the local section a : CP^ \ {oo} — )■ C^, given 
by [l,w] — ^ {^,w) and the linear map $ : — > which induces (p. Apply <& o ct 
to the point [l,w] to form 'i>((T([l, 77]) — (^,'^2) G and then apply the distance 
functions r^j to this configuration in to get the rij of eq. (18 1. That is, the 



functions of eq. ( 18 1 are r^j 0^0 a. Then 1 1^|| in the expression for V is the moment 
of inertia / = as given by ^ with the rij there being those given by eq. (18 1. 
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Figure 1 . Contour plot of the shape potential on the unit sphere 
s^+Sj+'^i = 1 in the equal mass case. There is a discrete symmetry 
of order twelve generated by the reflections in the sides of the 
indicated spherical triangle. 



Alternatively, we can view CP^ as S^/S^ and realize the S"^ by setting ||^|| = 1. 
Then the pij are r^j restricted to this S^, and then understood as invariant 
functions. 

Figure [T] shows a spherical contour plot of V for equal masses mi = m2 = rn^. 
The equator features the three binary collision singularities as well as three saddle 
points corresponding to the three collinear or Eulerian central configurations. The 
equilateral points at the north and south poles of the sphere are the Lagrangian 
central configurations which are minima of V . 

Figure [2] shows contour plots of the shape potential in stereographic coordinates 
(a;,y) for the equal mass case and for mi = l,m2 = 2, ma = 10. The unit disk 
in stereographic coordinates corresponds to the upper hemisphere in the sphere 
model. When the masses are not equal, the potential is not as symmetric, but due 
to the choice of coordinates, the binary collisions are still at the roots of unity and 
the Lagrangian central configuration (which is still the minimum of V) is at the 
origin. 

The variables in figures [T] and [2] are related by stereographic projection: 

_ 2x _ 2y _ 1 - - ?/2 

^ 1 + .T^ + j/^ ^ 1 + a;^ + y2 1 + + 

The following result about the behavior of the shape potential will be useful 
jl6j . Consider the potential in the upper hemisphere (the unit disk in stereographic 
coordinates). V{x,y) achieves its minimum at the origin. It turns out that V is 
strictly increasing along radial line segments from the origin to the equator. 
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Figure 2. Contour plot of the shape potential in stereographic 
coordinates {x,y). The unit disk corresponds to the upper hemi- 
sphere in the sphere model. On the left is the equal mass case as 
in figurejlj On the right, the masses are toi — 1, ni2 — 2, — 10. 



Proposition 4. ( Compare with lemma 4, section 6 of |16j .) For all positive masses, 
the shape potential V{x, y) satisfies 

xV^ + yVy = 0(a;, y)(l - - y^) 

where > with strict inequality if (x, y) ^ (0, 0). 

Proof. Write 

^ mim2 ^ rriim^ ^ TO2"^3 

Pl2 Pl3 P23 



with pij — rij/W^W, and r^, ||f|| expressed as functions of {x,y) using (18l and Q. 
Then a computation shows that 

xV^ + yVy = cl){x, y){l ~ x^ ~ y^) 

where 

(20) = -7 TTT-Tj (migi + 7712.92 + ^353) 

2(toi -I- 1712 + msjWQl 

and 

51 = ('^^a-^Dli--^) 52 - (r^3-^i2)(i-i) 93 ^ {rl3~rl3){-^ - -^). 

'12 '13 '12 '23 '13 '23 

Note that gi > with strict inequality except on the line where ri2 = ri^. Similar 
properties hold for 52 £^nd 53 and the proposition follows. QED 

2.3. Equations of Motion and Hill's Region. We can derive the equations of 
motion on Q = (0, 00) x CP^ by calculating the reduced Lagrangian L^ed in any 
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convenient coordinates and then writing out the resulting Euler-Lagrange equa- 



tions. We use the coordinates r,x,y as above with w — f?2/?7i ~ x + iy. (See (10 
^ and also eq. ([l8|, ([9]) and ([l9]). ) Then 



1-2, 1 n2,'-2i •2\ 3^i/X2 

+ ^K{x,y}r (x +y ) k=-^^-^. 



where ||^||^ as a function of x,y is obtained by plugging the expressions (18) into 
Lagrange's identity (|9|. Then 

Lred{r, X, y, r, x, y) = Kq + -V{x, y) 



r 



and so the Euler-Lagrange equations are 

1 



^2 V" + Kr{x^ + y 



2^ 



(21) [nr^xy = -V., + \K,r\x^+f) 

r I 



Conservation of energy gives 



K^-^-V{x,y) 
r 



Remark The expression K{x,y){dx'^ + dy"^) describes a spherically symmetric 
metric on the shape sphere. For example, when mi ~ TO2 ~ one computes that 
1^ — (xj^x'^j^y^Y which is the standard conformal factor for expressing the metric on 
the sphere of radius 1/2 in stereographic coordinates x, y. 

These equations describe the zero angular momentum three-body problem re- 
duced to 3 degrees of freedom by elimination of all the symmetries and separated 
into size and shape variables. An additional improvement is achieved by blowing 

3 

up the triple collision singularity at r = by introducing the time rescaling ' = r 2 ' 
and the variable v = r'/r [12_. The result is the following system of differential 
equations: 



(22) 



V =^v 



1 2 , / /2 , ,2 



k(x + y' ) — V 



(ny'y ^Vy^\nvy' + \Ky{x'^ +y'^). 



The energy conservation equation is 

(23) + \k{x''^ +y''^)-V{x,y) = -rh 

Note that {r = 0} is now an invariant set for the flow, called the triple collision 
manifold. Also, the differential equations for (v, x, y) are independent of r. Call the 
rescaled time variable s. Since the rescaling is such that t'{s) = r^ , behavior near 
triple collision that is fast with respect to the usual time may be slow with respect 
to s. Motion on the collision manifold could be said to occur in zero i-time. 

System (22) could be written as a first-order system in the variables (r, x,y,v, x' ,y'). 
Fixing an energy —h<0 defines a five-dimensional energy manifold 



Ph = {{r, X, y, V, x', y'):r>0, (23 1 holds}. 
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Figure 3. Half of the Hill's region for the equal raass three-body 
problem in stereographic coordinates {x, y, r) where the unit disk 
corresponds to the upper half of the shape sphere. In these coordi- 
nates, the configurations with collinear shapes are in the cylindrical 
over the unit circle. 

The projection of this manifold to configuration space is the Hill's region, Qh- Since 
the kinetic energy is non-negative the Hill's region is given by 

Qh^ {{r,x,y) -.Q <r <V{x,v)/h}. 

Figure [3] shows the part of the Hill's region over the unit disk in the (x,?/)-plane 
for the equal mass case. The resulting solid region has three boundary surfaces. 
The top boundary surface r = V{x,y)/h is part of the projection to configuration 
space of the zero-velocity surface, the bottom surface, r = 0, is contained in the 
projection of the triple collision manifold. The side walls are part of the vertical 
cylinder over the unit circle which represents collinear shapes. 

2.4. Visualizing the Syzygy Map. Figure |4] shows a different visualization of 
the same Hill's region. This time the shape is viewed as a point s — (si, S2, S3) on 
the unit sphere which is then scaled by the size variable r to form rs and plotted. 
The region of figure [3] corresponds to the upper half of the solid in the new figure. 
The collision manifold r = is collapsed the origin so that the bottom circle of 
the syzygy cylinder in figure [3] has been collapsed to a point. The collinear states 
within the Hill region now form the an unbounded "three-armed" planar region 
homeomorphic to the interior of the unit disk of figure [3] The syzygy map is the 
fiow-induced map from the top half of the boundary surface in figure |4] to the 
interior of this three-armed planar region. 

Figure [5] illustrates the behavior of the map in the equal mass case, using coor- 
dinates which compress the unbounded region Ch into a bounded one. The open 
upper hemisphere of the shape sphere has been identified with the domain of the 
syzygy map. The figure shows the numerically computed images of several lines of 
constant latitude and longitude on the shape sphere. The figure illustrates some 
of the claims of theorem [TJ Note that the image of the map seems to be strictly 
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Figure 4. The entire Hill's region for the equal mass three-body 
problem in coordinates r{si, S2, S3). The syzygy configurations Ch 
form the planar region, homeomorphic to a disk, dividing the Hill's 
region in half. The origin represents triple collision (which we count 
as a syzygy). The syzygy map maps from the upper half of the 
Hill region's boundary, to Ch- 



smaller than Ch and is apparently bounded by curves connecting the binary colli- 
sion points on the boundary. It does contain a neighborhood of the binary collision 
rays however. Also note that a line of high latitude, near the Lagrange homothetic 
initial condition (the North pole), maps to a small curve encircling the origin. Ap- 
parently there is a strong tendency, as yet unexplained, to reach syzygy near binary 
collision and near the boundary of the image. 

2.5. Flow on the collision manifold: linearization results. We will need some 
information about the flow on the triple collision manifold r = which we take from 



[T5] . Eq. ( 23 ) expresses the triple collision manifold as a two-sphere bundle over the 



shape sphere, with bundle projection (0, x, y, u, a;', y') {x, y). The vector field ( 22 ) 
restricted to the triple collision manifold r = flow has 10 critical points, coming 
in pairs one pair for each central configuration p. Two of these central 

configurations correspond to the equilateral triangle configurations of Lagrange and 
are located in our xy coordinates at the origin and at infinity, (x, y) = 00. We write 
L for the one at the origin. The other three central configurations are coUinear, 
were found by Euler, and are located on the unit circle in the xy plane, alternating 
between the three binary collision points. The two equilibria p± = {xo,yo,v±,0,0) 
for a given central configuration p ~ {xo,yo) are obtained by solving for v from 



the r = energy equation (231 to get v = v± = ±^y2V{xo, yo)- The positive 
square-root u+ corresponds to solutions 'exploding out" homothetically from that 
configuration, and the negative square-root v_ with w = u_ < corresponds to 
solutions collapsing into that central configuration. Associated to each central 
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Figure 5. Image of the syzygy map in the equal mass case. The 
open upper hemisphere of the shape sphere has been identified 
with the domain of the syzygy map. Several circles of constant lat- 
itude and a sector of arcs of constant longitude have been followed 
forward in time to the first syzygy and the resulting image points 
plotted. The coordinates on the image are arctan(r) (si, S2) so that 
the horizontal plane in figure [4] is compressed into the shaded open 
disk. The syzygy configurations Ch now form the region bounded 
by the heavy black curve. The image of the map, however, seems 
to be strictly smaller. 



configuration we also have the corresponding homothetic solution, which lives in 
r > and forms a heteroclinic connection connecting to p_ . 

We will also need information regarding the stable and unstable manifolds of the 
equilibria p±. This information can be found in |15j . Prop 3.3. 

Proposition 5. Each equilibrium p± — {O,xo,yo,v±,0,0) is linearly hyperbolic. 
Each has v± as an eigenvalue with corresponding eigenvector tangent to the homo- 
thetic solution and so transverse to the collision manifold. Its remaining 4 eigen- 
vectors are tangent to the collision manifold. 

The Lagrange point L_ has a 3 dimensional stable manifold transverse to the 
collision manifold and a 2-dimensional unstable manifold contained in the collision 
manifold. The linearized projection of the unstable eigenspace to the tangent space 
to the shape sphere is onto. At L+ the dimensions and properties of the stable and 
unstable manifolds are reversed. 

The Euler points i?_ each have a 2-dimensional stable manifold transverse to 
the collision manifold and contained in the collinear invariant submanifold and a 
3-dimensional unstable manifold contained in the collision manifold. At E^, the 
dimensions and properties of the stable and unstable manifolds are reversed. 
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Remarks. 1. The r > solutions in the stable manifold of i_ are solutions 
which limit to triple collision in forward time, tending asymptotically to the La- 
grange configuration in shape. The final arcs of the solutions of this set can be 
obtained by minimizing the Jacobi-Maupertuis length between points P and the 
triple collision point for P varying in some open set of the form U \ C where U 
is a neighborhoodof and C is the coUinear set. 

2. Because the unstable eigenspace for L_ projects linearly onto the tangent 
space to the shape sphere it follows that the unstable manifold of L_ cannot be 
contained in any shape sphere neighborhood a;^ + < e of L. 

3. The unstable manifold of an exploding Eulcr equilibrium lie entirely within 
the collinear space z — and real solutions lying in it form a two-dimensional set of 
curves. These curves can be obtained by minimizing the Jacobi-Maupertuis length 
among all collinear paths connecting points P and the triple collision point as P 
varies over collinear configurations in some neighborhood [/ H C of 0. 

4. The Sundman inequality implies that v' > everywhere on the triple collision 
manifold and that v{s) is strictly increasing except at the 10 equilibrium points. 
That is, V acts like a Liapanov function on the collision manifold. 



3. The Syzygy Map 

In this section the syzygy map taking brake initial conditions to their first syzygy 
will be studied. As mentioned in the introduction, it will be shown that every non- 
collinear, zero-velocity initial condition in phase space can be followed forward in 
time to its first syzygy. The goal is to study the continuity and image of the 
resulting mapping. 



3.1. Existence of Syzygies. In [THIIIT] Montgomery shows that every solution of 
the zero angular momentum three-body problem, except the Lagrange (equilateral) 
homothetic triple collision orbit, must have a syzygy in forward or backward time. 
In forward time, the only solutions which avoid syzygy are those which tend to 
Lagrangian triple collision. We now rederive this result, using our coordinates. 

The result will follow from a study of the differential equation governing the 
(signed) distance to syzygy in shape space. We take for the signed distance 



1-x' 



y 



the set of collinear shapes. From (21 1, one finds 



where x, y are the coordinates of section 2.3 Note the unit circle z = is precisely 



Kr^ 



Pi 



{xVx + yVy) - r^(i^ + y^)(K + ^{xk^ + yny))- 



A computation shows that 



K -|- ^(xKj. -|- yHy) — 



c(l — .T^ — 2/^) 



where 



3mim2»Ti3(miTO2 + mim^ H- m2m3) 
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Using this and proposition |4] gives 

. _ 2 

(24) ^ ~ K7^^' 

Pi = ~Fi{r,x,y,x,y)z 

where 

1 cr^(i^ + y^) 
Fi = -(^(.T,y) + . 

Note that Fi is a smooth function for r > and satisfies Fi > with equaUty only 
when x~y^x~y = 0. 



Using (24) one can construct a proof of Montgomery's resuh. First note that 
z = defines the syzygy set and z = pi = is an invariant set, namely the phase 
space of the coUinear three-body problem. Without loss of generality, consider an 
initial condition with shape (x, y) in the unit disk, i.e., the upper hemisphere in the 
shape sphere model. Our goal is to show that all such solutions reach z — 0. 



Remark. The key result of (TH] is a differential equation very similar to eq (24) 
for a variable which was also called z but which we will call zj now, in order to 
compare the two. The relation between the current z and this zj is 

z 

Zl = 



— z 

and can be derived from the expression z/ = ^j^^2jj^y2 for the height component of 
the stereographic projection map — > 5^ \ (0,0,-1). The important values of 
these functions on the shape sphere arc 

coUincar : z = 0; z/ = 

Lagrange, pos. oriented : z = 1; z/ = 1 

Lagrange, neg. oriented : z — +oo] zi ~ —1 

Proposition 6. Consider a solution with initial condition lying inside the punc- 
tured unit disc in the shape plane : < z(0) < ci < 1, and pointing outward (or 
at least not inward): z(0) < 0. Assume that the size of the configuration satisfies 
< r{t) < C2 for all time t > 0, and some positive constant ci. Then there is 
a constant To(ci,C2) > and a time to G [0,To] such that < z(t) < z(0) for 
t € [0,to) and z{to) = 0. 

Proof. Consider the projection of the solution to the (z,pi) plane. By hypothesis, 
the initial point (z(0),pi(0)) lies in the fourth quadrant of the plane. Since Fi > 0, 



(24 1 shows that z{t) decreases monotonically on any time interval [0,to] such that 



z(t) > so < z{t) < c\ holds. An upper bound for t^ will be now be found. 



Let a\ denote a clockwise angular variable in the (z,pi) plane. Then (24 1 gives 



2 1 2 

(z^ +|5^)ai = Fi{r,x,y,x,y)z'^ + — -pj > -<j){x,y)z^ H -pi 

Kr'' r Kr'^ 

2 nr^A n2 



On the set where < z < ci and < r < C2, the coefficients of z and pf of the 
right-hand side each have positive lower bounds. Hence there is a constant fc > 
such that cti > k holds on the interval [0,to]. It follows that to < 

It remains to show that the solution actually exists long enough to reach syzygy. 
It is well-known that the only singularities of the three-body problem are due to 
collisions. Double collisions are regularizable (and in any case, count as syzygies). 
Triple collision orbits are known to have shapes approaching either the Lagrangian 
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or Eulerian central configurations. The Lagrangian case is ruled out by tlie upper 
bound on z{t). Eulerian triple collisions can only occur for orbits in the invariant 
collinear manifold, so this case is also ruled out. QED 

The next result gives a uniform bound on time to syzygy for solutions far from 
triple collision. It is predicated on the well-known fact that if r is large and the 
energy is negative then configuration space is split up into three disjoint regions, 
one for each choice of binary pair, and within each region that binary pair moves 
approximately in a bound Keplerian motion. The approximate period of that mo- 
tion is obtained from knowledge of the two masses and the percentage of the total 
energy involved in the binary pair motion. The 'worst' case, i.e. longest period, is 
achieved by taking the pair to be that with greatest masses, in a parabolic escape 
to infinity so that all the energy —h is involved in their near Keplerian motion, 
and thus the kinetic energy of the escaping smallest mass is tending to zero. This 
limiting 'worst case' period is 

^ ^ 1 r 13/2 

* (2/i)3/2 ^mi-hm^-^ 

where the excluded mass mk is the smallest of the three. 
Here is a precise proposition. 

Proposition 7. Let be the constant above and let f3 be any positive constant less 
than 1. Then there is a (small) positive constant eo — £oW) such that all solutions 
with r(0) > energy —h < 0, and angular momentum have a syzygy within 
the time interval [0, ^r*]. Moreover, this syzygy occurs before r — |r(0), so that 
r{t) > at this syzygy. 

The final sentence of the proposition is added because if initial conditions are such 
that the far mass approaches the binary pair at a high speed then the perturbation 
conditions required in the proof will be violated quickly: r(t) will become 0(1) in a 
short time, well before the required syzygy time ^t*. In this case the approximate 
Keplerian frequency of the bound pair is also accordingly high, guaranteeing a 
syzygy well before perturbation estimates break down and well before the required 
syzygy time. 

Proof. The proof is perturbation theoretic and divides into two parts. In the first 
part we derive the equations of motions in a coordinate system quite similar to the 
one which Robinson used to compactify the infinity corresponding to r — > oo at 
constant h. In the second part we use these equations to derive the result. 

Part 1. Deriving the equations in the new variables. We will use coordinates 
adapted to studying the dynamics near infinity which are a variation on those intro- 
duced by McGehee ([13 ) and then modified for the planar three-body problem by 
Easton, McGehee and Robinson ( [SJ [3 [20] ) . Going back to the Jacobi coordinates 
Ci,6 set 

thus defining coordinates {p,u,9) e x C x S^. The variables p,u coordinatize 
shape space while 9 coordinatizes the overall rotation in inertial space. Then 



2 I |2 I 2 
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The reduced kinetic energy (i.e. metric on shape space) is given by 
Here 

u Ail = U1U2 — U2U1 = Im{uu). 

(Kq can be computed two ways: either plug the expressions for the in terms of 
p,u,9,... into the expression for the kinetic energy and then minimize over 0, or 
plug these same expressions into the reduced metric expression.) Next we make a 
Levi-Civita transformation by setting u = which gives 

Ko = 2f,,\zf\z\' + _ ^^(^ ^ ^)2_ 

We want to introduce the conjugate momenta and take a Hamiltonian approach. 
We can write 

Ko = 2ij.i\zf{A{p, z))z,z) + -ii2f)^ 
where A is the symmetric matrix A = I — B with 

|2 



Z2 -Z1Z2 
-Z1Z2 zl 

Here Z1.Z2 denote the real and imaginary parts of z, not new complex variables. 
The conjugate momenta are 

r] = 'iiii\z\^ Az y ^ H2p. 

If we write 

A-^ = I + ,i^f{p,z) 

we find 

For later use, note that p,if = B + B"^ + . . . is a positive semi-definite symmetric 

matrix and / = 0(|z|'*/r^). 

The negative of the potential energy is 

mim2 mp2 , , x 
U{p,z) = — ^ + +g{p,z) 

where the "coupling term" g is 

9(P^^) = TT^r 211 + n 211 =Oi zf p^) 

11/0 + t^2Z II \\P - VlZ^W P VI I / / 

For r > tq sufficiently large the Hill region breaks up into three disjoint regions. 
We are interested for now in the region centered on the 12 binary collision ray. In 
this case, for r > ro sufficiently large one finds that |2;| is bounded by a constant 
depending only on the masses and h and the choice of tq. This bound on |2;| tends to 
a nonzero constant as ro — *• 00. It then follows from the identity = p,\p'^ + p,2\z\^ 
that r = y/p2p + 0(1) for r large. 

Next we complete the regularization of the binary collision by means of the 
the time rescaling ^ = l^^l^^- Using the Poincare trick, the rescaled solutions 
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(25) 



with energy —h become the zero-energy solutions of the Hamiltonian system with 
Hamiltonian function, 

H=\z\Ko-U + h) 

= + \^y' + lifiP, zH V) - m,m, - - g{p, .)|.|^ + h\z\\ 

0/J.i Z/i2 O p 

Computing Hamihon's equations, making the additional substitution 

X = 1/p 

to move infinity to a; = 0, and writing ' for d/dr gives the differential equations: 

I^P 2 

X = — ' X y 

M2 

y' = -mfi2\zfx^ + \z\'^gp - -{fpV,v) 

z' = -j—V + \fr} 
4/ii 4 

1 

r?' = -2(^+ ^ - miJ,2X - g)z - - {fzV,v) + \z\'^9z 

where the subscripts on / and g denote partial derivatives. These satisfy the 
bounds: 

f = 0{xM'') h = 0{xM'') fp = 0{xM^) 

g = 0{xM^) g, = 0{xM^) gp = 0{xM% 

The energy equation is H = 0. 

Infinity has become a; = 0, an invariant manifold. At infinity we have x' =y' = Q 
while other variables satisfy 

, 1 

z = -. — n 

ri = -2{h+^)z. 

Since y is constant this is the equation of a two-dimensional harmonic oscillator. 

Part 2. Analysis. Observe that the configuration is in syzygy if and only if 
the variable u is real . Since u = we have syzygy at time t if and only if z{i) 
intersects either the real axis or the imaginary axis. If ^;'s dynamics were exactly 
that of a harmonic oscillator then it would intersect one or the other axis (typically 
both) twice per period. We argue that z is sufficiently close to an oscillator that 
these intersections persist. 

We begin by establishing uniform bounds on z and r\ valid for all x sufficiently 
small. These come from the energy. We rearrange the expression for energy into 
the form 

1 IzP 1 

7, — + 7. — + 7;{.f{P,z)v,v) + {h- mp2X - g{x,z))\z\'^ = mima. 

opi Zp2 o 

Choose Co so that x < cq implies \mp2X+g\ < h/2. By the positive semi-definiteness 
of / we have 

^\r]\^ + {h/2)\z\^ < mim2 
which gives our uniform bounds on z,r]. 
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It now follows from equations ( 25 ) that for all e < eg and x < e we have 

\x'\ < Cx^\y\ 



\y 

(26) , 1 

z = — 

Am 

Tj' = 2Hi2Z + 0{e 



< Cx^ 



where H12 = ~h — + mfj,x + g, a quantity which represents the energy of the 
binary formed by masses mi, 7712. Here C is a constant depending only on the 
masses, the energy and Eq. 

Fix an initial condition xo,yo, zq^tiq with xq < e/2. Write {x{t),y{t), z(t),rj(t)) 
for the corresponding solution. For C a positive constant, write TZ = TZc,e for the 
rectangle < a; < e, jy — ?/o| < Ce in the xy plane. We will show there exists Ci 
depending only on the masses, h, on the constant C and j/o such that the projection 
x{t),y(t) of our solution lies in TZc,e for all times t with \t\ < 1/Cie. Suppose that 
the projection of our solution leaves the rectangle TZ in some time t. If it first leaves 
through the y-side, then we have Ce = \y{t) — yo\ = \J ydt\ < Ce^ J dt = Ce^t, 
asserting that 1/e < t. Thus it takes at least a time 1/e to escape out the y-side. 
To analyze escape through the x-side at x = e we enlist Gronwall. Let 6*2(1/0) = 
C'dyol + Ce). Compare x{t) to the solution x to i ~ G2X^ sharing initial condition 
with x{t), so that a;(0) = xq. The exact solution is i = xq/{\ — XQC2t). Gronwall 
asserts x{t) < x{t) as long as x{t),y{t) remain in the rectangle (so that the estimates 
(26 1 are valid). But x{t) < e ioi t < l/eC2- Consequently it takes our projected 
solution at least time t — l/eC2 to escape out of the x-side, and thus {x{t),y{t)) 
lies within the rectangle for time \t\ < 1/eCi, with 1/C\ = min{l, l/C2(yo)}- 

We now analyze the oscillatory part of equations (25). We have 

2 2 
H12 ^ -h- h m^x + g = -h - h 0{e). 

2^2 2^2 

and let H'^2 — ~ ^Vo- Then as long as (x, y) £ TZc.e we have the bound 

\H',2~H,2{x,y,z,7j)\^0{e). 

holds. Thus the difference between the vector field defining our equations and the 
"frozen oscillator" approximating equations 

, 1 

z — 77 

(27) 4/^1 

V' = 2H°2Z 



which we get by throwing out the 0(e ) error terms in (26) and replacing H12 by 
i??2 is Oie). 

Now the period of the frozen oscillator is Tq = 2n^/2jl^/ y/\Hf^\. On the other 
hand we remain in TZc.e for a-t least time T-jz = 1/eCi. Both of these bounds depend 
on yo but we have 

Tq ^ 27r^/27I7max{l, C2{yo)} ^ ^ 
Tn ' Vh + ym ~ 
where C3 does not depend on y^. Moreover we have a uniform upper bound 
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Hence for eo > sufficiently small and xq < e/2, e < eg, a solution of (26) 
remains in TZc.e for at least time Tmax a-nd the difference between the component 
z{t) of our solution and the corresponding solution z to the linear frozen oscillator 



equation (27) is of order e in the C^'^-norm: because we also get the 0(e) bound 
on 77 = A^iz' . Now, any solution z to the frozen oscillator crosses either the real 
or imaginary axis, transversally, indeed at an angle of 45 degrees or more, once per 
half-period. (The worst case scenario is when the oscillator is constrained to a line 
segment). Thus the same can be said of the real solution z{t) for eg small enough: 
it crosses either the real or imaginary axis at least once per half period. 

Finally, the time bounds involving of the proposition were stated in the New- 
tonian time. To see these bounds use the inverse Levi-Civita transformation, and 
note that a half-period of the Levi-Civita harmonic oscillator corresponds to a full 
period of an approximate Kepler problem with energy i?"2- This approximation 
becomes better and better as e — 0. Moreover, the Kepler period decreases monon- 
tonically with increasing absolute value of the Kepler energy so that the maximum 
period corresponds to the infimum of |-ffi2| and this is achieved by setting ?/ = 
(parabolic escape), in which case — ~h and we get the claimed value of r* with 
/? = 1. (Letting /3 — > 1 corresponds to letting eg 0. ) QED 

These two propositions. Proposition [7] and Proposition [6| lead to Montgomery's 
result: 

Proposition 8. Consider an orbit with initial conditions satisfying < z(0) < 
1 and r(0) > 0. Either the orbit ends in Lagrangian triple collision with z{t) 
increasing monotonically to 1, or else there is some time > such that z{ta) — 0. 

Proof. Consider a solution with no syzygies in forward time. By proposition [7] 
there must be an upper bound on the size: r(t) < C2 for some C2 > 0. Next suppose 
z(0) = 1, i.e., the initial shape is equilateral which means {x{0),y{0)) = (0,0). If 
(i(0),y(0)) = (0,0) the orbit is the Lagrange homothetic orbit which is in accord 
with the proposition. If {x{0),y{0)) 7^ (0,0) then for all sufficiently small positive 
times ti one has < z{ti) < 1 and z{ti) < 0. Then proposition [6] shows that there 
will be a syzygy. It remains to consider orbits such that < z(0) < 1. 

Suppose < z(0) < 1 and that the orbit has no syzygies in forward time. It 
cannot be the Lagrange homothetic orbit and, by the first part of the proof, it 
can never reach z{t) = 1. Proposition [6] then shows that z{t) > holds as long 
as the orbit continues to exist, so z{t) is strictly monotonically increasing. By 
proposition [7] there is a uniform bound r{t) < C2 for some 02- There cannot be a 
bound of the form z(t) < Ci < 1, otherwise a lower bound di > A: > as in the 
proof of proposition [6] would apply. It follows that z{t) — > 1. 

To show that the orbit tends to triple collision, it is convenient to switch to the 



blown-up coordinates and rescaled time of equations (22). Since the part of the 
blown-up energy manifold with z > z{0) > is compact, the w-limit set must be a 
non-empty, compact, invariant subset of {z = 1}, i.e., {{x,y) = (0,0)}. The only 
invariant subsets are the Lagrange homothetic orbit and the Lagrange restpoints 
L±. Since W^{L+) C {r = 0}, the orbit must converge to the restpoint L_ in the 
collision manfold as s — > 00. In the original timescale, we have a triple collision 
after a finite time. QED 

3.2. Existence of Syzygies for Orbits in the Collision Manifold. In the last 
subsection, it was shown that zero angular momentum orbits with r > and not 
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tending monotonically to Lagrange triple collision have a syzygy in forward time. 
This is the basic result underlying the existence and continuity of the syzygy map. 
However, to study the behavior of the map near Lagrange triple collision orbits we 
need to extend the result to orbits in the collision manifold {r = 0}. This entails 
using the equations (22 1. 

Setting z — 1 — x'"' — as before one finds 

, 2 

Z = -P2 



K 

P2 = -{xVx + yVy) - {x'"^ + y''^){K + \{XK,J, + yKy)) + \VP2. 

As in the last subsection, this can be written 

(28) ^ l^' 

p'2 = -F2{x,y,x',y')z - \vp2 

where 

F2^4>{x,y)+c{x'^ + y'^). 

F2 is a smooth function satisfying F2 > with equality only when x = y = x' = 
y' = 0. Recall that ' denotes differentiation with respect to a rescaled time variable, 
s. 

To see which triple collision orbits reach syzygy, introduce a clockwise angle a2 
in the {z,p2)-'p\&'ae. Then 

(29) (z^ + p^a'2 = F2Z^ + -pI + \vzp2. 

The cross term in this quadratic form introduces complications and may prevent 
certain orbits with r = from reaching syzygy. 

To begin the analysis, consider a case where the cross term is small compared 
to the other terms, namely orbits with bounded size near binary collision. 

Proposition 9. Let C2 > 0. There is a neighborhood U of the three binary collision 
shapes and a time Sq{U,C2) such that if an orbit has shape {x{s),y{s)) S U and 
size r{s) < C2 for s G [0,6*0] then there is a (rescaled) time sq G [0,5*0] such that 
z(so) =0. 

Proof. Let S, t be the determinant and trace of the symmetric matrix of the qua- 
dratic form 

Iv 2k^i 

Iff (5, r > then the smallest eigenvalue of this matrix is greater than 2S/t, which 
yields the lower bound 

, 2(5 „ 
ao > — > 0. 

^ ~ T 



From ( 18 ) it follows that near binary collision, one of the three shape variables 
w while the other two satisfy r^^, r^-^ ~ \Pi. Consider, for example, the binary 
coUision r\2 = at (x,?/) = (1,0). One has F2 > 4'{x,y) and (j){x,y) « C/rfj for 
some constant C > 0. The matrix entry is bounded and v can be estimates 
using the energy relation ( 23 1 

< 2{V{x, y) - rh) = 2mim2/ri2 + 0(1) 
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near {x,y) = (1,0), r < C2. Using these estimates in the trace and determinant 
gives the asymptotic estimate 

Similar analysis near the other binary collision points yields a neighborhood U in 
which there is a positive lower bound for This forces a syzygy (z = 0) in a 
bounded rescaled time, as required. QED 

This result, together with some well-known properties of the flow on the triple 
collision manifold leads to a characterization of possible triple collision orbits with 
no syzygy in forward time. 

Proposition 10. Consider an orbit with initial condition r{0) — and < z(0) < 
1. Either the orbit tends asymptotically to one of the restpoints on the collision 
manifold or there is a (rescaled) time sq > such that z(so) — 0. 



Proof. Recall from section |2.5| that the flow on the triple collision manifold is 
gradient-like with respect to the variable v, i.e., v{s) is strictly increasing except 
at the restpoints. It follows that every solution which is not in the stable manifold 



of one of the restpoints satifies v{s) — )■ 00. In this case, the energy equation (23) 
shows that the shape potential V{x{s), y{s)) 00 so the shape must approach one 
of the binary collision shapes. For such orbits, proposition |9] gives a syzygy in a 
bounded rescaled time and the proposition follows. QED 

We will also need a result analogous to proposition |6j Consider an orbit in 
{r — 0} with initial conditions satifying < z(0) < ci < 1 and z'{0) < 0. It 
will be shown that, under certain assumptions on the masses, every such orbit has 
z{sq) = at some time sq > 0. 

On any time interval (0, sq] such that z{s) > one has z'{s) < and hence 
z{s) is monotonically decreasing. This follows from the "convexity" condition that 



z" < whenever < z < 1 and z' = which is easily verified from (28). If the 
orbit does not reach syzygy in forward time, then < z{s) < ci for all s > 0. 
It is certainly not in the stable manifold of one of the Lagrangian restpoints at 



z = 1 so, by proposition 10 it must be in the stable manifold of one of the Eulerian 
(collinear) restpoints at z = 0. 

It will now be shown that, for most choices of the masses, even orbits in these 
stable manifolds have syzygies. Let Cj = (xj,yj),j = 1,2,3 denote the collinear 
central configuration with mass ruj between the other two masses on the line. The 
two corresponding restpoints on the collision manifold are Ej- and Ej^ with co- 
ordinates {r,x,y,v,x' ,y') = (0, Xj,yj, ±y^2V{xj ,yj), 0, 0). The two-dimensional 
manifolds ^{Ej-) and W'^{Ej+) are contained in the collinear invariant subman- 
ifold. In particular, an orbit with z(0) > 1 cannot converge to Ej- in forward 
time. On the other hand, W'^{Ej+) is three-dimensional and its intersection with 
the collision manifold is two-dimensional. These are the orbits which might not 
reach syzygy. 

Remark. Conjecture [l] asserts that all orbits reach syzygy with v < 0. So, if the 
conjecture is valid then convergence to Ej^ would be impossible and the discussion 
to follow, and the restriction on the masses in the theorem, would be unnecessary. 

For most choices of the masses, the two stable eigenvalues of Ej^ with eigenspaces 
tangent to the collision manifold are non-real. See [TS] . In this case we will say that 
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Ej^ is spiraling. The spiraling case is more common: real eigenvalues occur only 
when the mass ruj is much larger than the other two masses. For all masses, at 
least two of the three Eulerian restpoints Ej^ are spiraling and for a large open set 
of masses where no one mass dominates, all of the Eulerian restpoints are spiraling. 

Proposition 11. Let Ej^ he a spiraling Eulerian restpoint. Then there is a neigh- 
borhood U of Ej-f-, and a time So{U) such that any non-collinear orbit with initial 
condition in the local stable manifold W^{Ej^) has a syzygy in every time interval 
of length at least Sq . 

Proof. Introduce local coordinates in the energy manifold near Ejj^ of the form 
{r,a,b, z,p2) where (a, 6) are local coodinates in the collinear collision manifold 
(the intersection of the invariant collinear manifold with {r = 0}) and (2,^2) are 



the variables of ( 28 1 . The invariance of the collinear manifold z ^ p2 = implies 



that the linearized differential equations for (2,^2) take the form 



z 


1 


'a (3 




z 






7 S 







where the eigenvalues of the matrix are the non-real eigenvalues at the restpoint. 
The spiraling assumption implies that the angle 012 in the (2;,p2)-plane satisifies 
a'2 > k > Q for some constant k. So for the full nonlinear equations, one has 
a2 > fc/2 > in some neighborhood U . Since non-collinear orbits in the local 
stable manifold remain in U and have nonzero projections to the (z,p2)-plane, the 
proposition holds with 6*0 = 2-k jk. QED 

A mass vector (rrii, m2, TO3) will said to satisfy the spiraling assumption if all of 



the Eulerian restpoints are spiraling. The next result follows from propositions 10 
andini 

Proposition 12. Consider an orbit in {r — 0} with initial conditions satifying 
< z(0) < Ci < 1 and z'(0) < 0. Suppose that the masses satisfy the spiraling 
assumption. Then there is a (rescaled) timetime sq > such that z{so) = 0. 
Moreover z{s) is monotonically decreasing on [0, sq] 

3.3. Continuity of the Syzygy Map. In this section, we will prove the statement 
in Theorem [1] about continuity of the syzygy map and its extension to the Lagrange 
homothetic orbit. 

We begin by viewing the syzygy map in blow-up coordinates {r,x,y,v,x' ,y'). 
Recall the notations Ph for the energy manifold and Qh for the Hill's region. Points 
on the boundary surface dQh of the Hill's region can be uniquely lifted to zero- 
velocity (brake) initial conditions in Ph- Let dQ^j^ be the subset of the boundary 
with shapes in the open unit disk (which corresponds to the open upper hemisphere 
in the shape sphere model). This is the upper boundary surface in figure [Sj It will 
be the domain of the syzygy map. 



To describe the range, let Ch = {(r, x, y,v,x' ,y') : r > OjX^+y"^ — 1, ( 23 ) holds} 
be the subset of the energy manifold whose shapes are collinear and let Ch C Qh 
be its projection to the Hill's region, i.e., the set of syzygy configurations having 
allowable energies. Ch is the cylindrical surface over the unit circle in figure [s] (the 
unit circle in {r = 0} is the blow-up of the puncture at the origin) . The first version 
of the syzygy map will be a map from part of to Ch. 

Recall that the origin {x, y) = (0, 0) represents the equilateral shape. The cor- 
responding point po £ (^Qh ^^^^^ ^ brake initial condition in Ph which is on the 
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Lagrange homothetic orbit. This orbit converges to the Lagrange restpoint without 
reaching syzygy. It turns out that pq is the only point of dQ'^ which does not reach 
syzygy and this leads to the first version of the syzygy map. 

Proposition 13. Every point of dQ^ \ po determines a brake orbit which has a 
syzygy in forward time. The map F : dQ^ \ Po ^ determined by following these 
orbits to their first intersections with and then projecting to Ch is continuous. 

Proof. The surface dQh intersects the line x ~ y = Q transversely at po. So every 
point of dQ^ \ po satisfies < z < 1 where z = 1 — — as before. Also r > 
on the whole surface dQh. 

Consider the brake initial condition corresponding to such a point. Since all 
the velocities vanish one has i(0) = 0. As in the first paragraph of the proof of 
proposition [6j one finds that z{t) is decreasing as long as z(t) > 0. In particular, 
z{t) does not monotonically increase toward 1. It follows from proposition [s] that 
there is a time to > Q such that z{tQ) — 0. In other words, the forward orbit reaches 

To see that the flow-defined map F : dQ^ \ Po ^ Ch is continuous, note that 
if the first syzygy is not a binary collision then z{to) < 0. This follows since z{t) 
is decreasing and since the orbit does not lie in the invariant coUinear manifold 
with z = i = 0. So the orbit meets Ch transversely. After regularization of 
double collisions, even orbits whose first syzygy occurs at binary collision can be 
seen as meeting Ch transversely as in the proof of proposition [A It follows from 
transversality that F is continuous. Composing with the projection to the Hill's 
region shows that F is also continuous. QED 

To get a continuous extension of F to po: we have to collapse the triple collision 
manifold back to a point. One way to do this is to replace the coordinates (r, x, y) 
with {X, Y, Z) — r(si, S2, ss) where s, are given by inverse stereographic projection. 
The Hill's region is shown in figure 4] In this figure, dQ'^ is the open upper half 
of the boundary surface and Ch is the planar surface inside (minus the origin). In 
this model, triple collision has been collapsed to the origin (X, Y, Z) = (0, 0, 0). It 
is natural to extend the syzygy map by mapping the Lagrange homothetic point pq 
to the triple collision point, i.e., by setting F{pq) = (0,0,0). The extension maps 
into Ch^ChUO. 

Theorem 5. If the masses satisfy the spiraling assumption then the extended syzygy 
map F : dQ^ — !■ Ch is continuous. Moreover, it has topological degree one near po . 

Proof. To prove continuity at Pq it suffices to show that points in dQ^ near pq have 
their first syzygies near r = 0. The proof will use the blown-up coordinates (r, x, y) 
and the rescaled time variable s. 

Let L_ denote the Lagrange restpoint on the collision manifold to which the 
orbit of Pq converges. Let iy^+(L_), VF"+(L_) be the local stable and unstable 
manifolds of L_ where local means that the orbits converge to L_ while remaining 



in {z > 0}. It follows from proposition 12 that z{s) is monotonically increasing 



to 1 along orbits in W'^~^{L-) D {r = 0}. A similar argument with time reversed 
shows that z{s) monotonically decreases from 1 along orbits in T/F"+(L_). These 
last orbits lie entirely in {r = 0}) by proposition [5] 

The unstable manifold W"'^{L^) is two-dimensional and its projection to the 
{x,y) plane is a local diffeomorphism near L_. Let be a small disk around L_ 
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in It follows from proposition 12 that all of the points in 13" \ L„ can 



be followed forward to meet Ch transversely. Moreover the monotonic decrease of 
z{s) implies that if this flow-induced mapping is composed with the projection to 
Ch and then to the unit circle, the resulting map from the punctured disk to the 
circle will have degree one. 

Let 7"" = dD'^ be the boundary curve of such an unstable disk. Since the unstable 
manifold is contained in the collision manifold, the first-syzygy map takes 7" into 
C*/! n {r = 0}. Transversality implies that the first syzygy map is defined and 
continuous near 7". Hence, given any e > there is a neighborhood U of 7u such 
that initial conditions in U have their first syzgies with r < e. Standard analysis of 
the hyperbolic restpoint L_ then shows that any orbit passing sufficiently close to 
L_ will exit a neighborhood of L_ through the neigborhood U of 7" and therefore 
reaches its first syzygy with r < e. We have established the continuity of the syzygy 
map at po- 

Now consider a small disk, D, around po in ^Q'^. It has already been shown 
that every point in D \ pq can be followed forward under the fiow to meet Ch 
transversely with z{s) decreasing monotonically. If D is sufficiently small it will 
follow the Lagrange homothetic orbit close enough to for the argument of the 
previous paragraph to apply. In other words the syzygy map takes D \ pq into 
{r < e} which implies continuity of the extension at pq. As before, the monotonic 
decrease of z{s) implies that the map from the punctured disk to the unit circle has 
degree one. Hence the map with triple collision collapsed to the origin has local 
degree one. QED 

This completes the proof of the continuity statements in Theorem [T] and also 
shows that the extended map covers a neighborhood of the triple collision point. 
See figure [5] for a picture illustrating this theorem in the equal mass case. 

4. VARIATIONAL METHODS: EXISTENCE AND REGULARITY OF MINIMIZERS 

In this section, we prove theorems [3] |4] and lemmas [T] by applying the direct 
method of the Calculus of Variations to the Jacobi-Maupertuis [JM] action: 

(30) A/m(7)= / ^/KoVAU-h) dt. 

Here 7 : [to,ti] — > Qh is a curve. Curves which minimize Ajm within some class of 
curves will be called "JM minimizers" . 

The integrand of the functional Ajm is homogeneous of degree one in velocities 
and so is independent of how 7 is parameterized. The natural domain of definition 
of the functional is the space of rectifiable curves in Qh- We recall some notions 
about Frechet rectifiable curves. See [8] for more details. Kinetic energy Kq induces 
a complete Riemannian metric, denoted 2Kq on the manifold with boundary Qh- 
We denote by dp the associated Riemannian distance. The Frechet distance between 
two continuous curves 7 : [to,ti] — > Qh and 7' : [io,i'i] Qh is defined to be 

p(7,7') = inf sup do(7(t),7'(Mi))), 
te[to,ti] 

where the infimum is taken over all orientation preserving homeomorphisms h : 
[to,ti] — >■ [io,ii]. We consider two curves equivalent if the Frechet distance between 
them is zero. An equivalence class of such curves can be seen as an unparametrized 
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curve. The Frechet distance is a complete metric on the set of unparametrized 
curves. The length £(7) of a parametrized curve 7 : [tQ, ti] Qii can be defined to 
be the supremum of the lengths of broken geodesies associated to subdivisions of 
[to,ii]- This length equals the usual 2Ko length when the curve is C^. Equivalent 
curves have equal lengths. A curve is called rectifiable if its length is finite. For a 



generic parametrization of a rectifiable curve, the integral (30) must be interpeted 
as a Weierstrass integral (see [5|). This integral is independent of parameterization. 

We will use the following two facts in our proof of Theorem [3] If -ftT C Q/i is 
compact, then the set of curves in Qh which intersect K and have length bounded 
by a given constant forms a compact set in the Frechet topology. (This fact is 
a theorem attributed to Hilbert.) The second fact asserts that both the length 
functional C and the JM action functional Ajm are lower semicontinuous in the 
Frechet topology. (Again, see 8 ). 

Since the Jacobi-Maupertuis action degenerates on the Hill boundary we will 
need to construct a tubular neighborhood of the Hill boundary within which we 
can characterize JM-minimizers to the Hill boundary. Our construction follows 
[221. 



Proposition 14. There exists e > 0, a neighborhood Sf^ of dQh (in Qh), o,n (an- 
alytic diffeomorphism $ : dQh x [0, e] — >■ Sf^ (satisfying ^{dQii,0) = dQh) and 
a strictly positive constant M, such that if 5 ^ (Oj^); ^'^'^ ^ £ dQh, the curve 
y 1—^ *&(x, y), y (z [0, S] is a reparametrization of an arc of the brake solution start- 
ing from X, and its length is smaller or equal to MS. If q ^ $(a;,(5) then for every 
rectifiable curve 7 in Qh joining q to dQh we have 

with equality if and only if ^ is obtained by pasting ($(a;, 2/))j^gjQ to any arc con- 
tained in dQh- Moreover, there exists a strictly positive constant a, independent 
from 5, such that U > h + a5 on Qh \ Sh, where we term Sh — ^{dQh, [0, S]). 

We postpone the proof of this Proposition, and instead use it now to prove 
Theorem [3] and Lemma [l] We will say that Sf^ is a Seifert tubular neighborhood of 
the Hill boundary, and that ^{dQh, e) is the inner boundary of Sf^. We prove now 
the first part of Theorem |3] 

Proposition 15. Given a point go in the interior of the Hill region Qh, there exist 
a JM action minimizer among rectifiable curves joining qo to the Hill boundary 
dQh- 



Proof. Let S'^ be the Seifert tubular neighborhood given by Proposition 14 If 
qo G S"^ then the unique brake orbit joining qg to the Hill boundary is a JM- 
minimizer. 

If qo ^ 5^, let 7„ be a JM minimizing sequence of rectifiable curve joining qo 
to dQh- Without loss of generality we can assume all these curves are Lipschitz 
and defined on the unit interval [0, 1]. Let i„ G (0, 1) be the first time such that 
7„(f) touches the inner boundary $(9Q/i,e) of S^. Let us replace 7 |[t„,i] by the 



unique brake solution joining ^n(tn) to dQh- By Proposition 14 this modification 



decreases the JM action of 7„, so our sequence is still a minimizing one. Let C be 



an upper bound for the numbers AjMi'^n)- Applying Proposition 14 again we get 
(31) C > AjMhn) > e'/' + C{^n |[0,t„] 
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therefore 

C - e3/2 

(32) £(7„) < —— + Me, 



for all n, proving that the lengths C{jn) are bounded. Since 7n(0) = qo for all 
n, Hubert's theorem discussed above applies: the sequence of curves 7„ is rela- 
tively compact in the Frechet topology. Therefore a subsequence of this sequence 
converges to a rectifiable curve 7 in Qh- This 7 is a JM minimizer by the lower 
semicontinuity of the JM action. QED 

Let d J m{<1o ill) denote the infinium of the JM-action among rectifiable curves 
in Qh joining go to gi, and let djAiilo, dQh) denote the minimum of the JA/-action 
among rectifiable curves in Q/j joining go to the Hill boundary. 

Lemma [1] the JM Marchal lemma, proof. 

Proof. First we establish the existence of a minimizer. 

If djMiqo,qi) > djMiqo,dQh) + djMiqi,dQh), then take a curve realizing the 
minimum in djM{qoidQh), another curve realizing the minimum in djM{qi,dQh) 
and join these two curves by any curve lying on the Hill boundary and connect- 
ing the endpoints of these two. In this way we get a rectifiable curve 7 join- 
ing go to gi and (possibly) spending some time on the Hill boundary. Since the 
JM action of any curve on the Hill boundary is zero, we have that Ajm{i) = 
djM{qo,dQh) + djM{qi,dQh)- Hence 7 is our minimizer. (As a bonus we have 
shown that djM{qo,qi) = djMiqo,dQh) + djM{qi,dQh) in this case.) 

If djM{qo,qi) < djM{qo,dQh) + djM{qi,dQh), let 7„ be a JAf- minimizing se- 
quence of rectifiable curves in Qh joining go to gi. We show now that there exists 
e > such that for n is sufficiently great, the 7„ do not intersect the Seifert 
tubular neighborhood S*^. Assume, for the sake of contradiction that there in 
fact exists a decreasing sequence of strictly positive real numbers e„ such that 
e„ — )■ as n -l-oo and a subsequence of the minimizing sequence, still denoted 
(7„ : [0, 1] — ?► Qh)n, and a sequence of times t„ such that for every n, 7„(t„) is 
in S'j^ . Let 7„ be the curve defined in the following way. Follow 7„ from go to 
ln{tn)- Join 7„(t„) to dQh by the unique brake solution in S^^ with one end point 
7„(t„). Return along the same brake solution, to 7n(^n). Then continue along 7„ 
from 7„(i„) to gi. The curve 7„ joins go to gi is rectifiable and touches dQh, hence 



AjMiln) > djMiqo,dQh) + djM{qi,dQh)- By proposition 14 we have 



djMiqo,dQh)+djMiqi,dQh) - 2ef/2 < AjMijn) ~ 2e3/2 < AjMiln)- 

Taking the limit n -> -l-oo we get djM{qo,dQh) + djM{qi,dQh) < rfjA/(<?o, gi), a 
contradiction. We have established our desired e > with the property that none 
of the curves 7„ {n sufficiently large) intersect S*^. 

Let C > be an upper bound of AjmIju) and a — a{e) the positive constant 



given by Proposition 14 A computation similar to (31) and (32) gives 



Hubert's Theorem again yields a curve 7 in Qh joining go to gi which the a subse- 
quence of the 7„ converges to in the Frechet topology. By lower semicontinuity of 
AjM, the curve 7 is a JM-minimizer. 
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Finally, we show that any minimizer 7 : [Iq, ti] — )• Qh is collision-free and that ev- 
ery subarc of 7, upon being reparametrized, is a true solution to Newton's equation. 
The classical Lagrangian action on the reduced space Q is defined by 

Aed(o-) / Lred{cr{t),&{t))dt 
Jto 

where a : [to,ti] ^ Q is any absolutely continuous path and the reduced Lagrangian 



Lred is defined in (14). lit € [to, ti] is such that 7(i) G dQh, obviously 7(i) is not a 
collision, hence it suffices to prove the statement for subarcs 7 \ [a,b] , [oj C [ioj^i] 
of 7 which never touch the Hill boundary. Let [a, b] be such an interval. By the 
inequality + > 2AB applied to the factors A — \/Ko, B — \/\J — h of the 
integrand for Aj]\j we see that 

V2AjM (7 I [aM] ) < Aredil \ [aM] ) ~ h{b ~ a) 

with equality if and only if the energy Kq — U evaluated on 7 is equal to —h 
for almost all t £ [a,b]. We may assume, without loss of generality, that 7 is 
parametrized by its 2iiro-arclength parameter s so that in the integrand Kq{s) = 1/2 
for almost every s e [a, 6]. Since 7 is a JM-minimizer, we know that 7(a) 7^ "fib), 
that the action Ajm{j \ [a,b] ) is finite, and that the set of s € [a, b] such that 7(5) 
is a collision is a closed set of zero Lebesgue measure. Let a : [a, /3] — )■ Qh be the 
reparametrization of 7 \ [a.b] defined by a{t) = j{s{t)), where s = s{t) is the inverse 
function of 

t = t{s)= [ '^^ se[a,6]. 

J a Vmiliv)) - h) 

The function t = t{s) is C^, strictly increasing and ^(s) > on an open set with 
full measure, hence the inverse function s = s{t) is absolutely continuous, and since 
s I— > 7(5) is Lipschitz, t i— > a{t) = ^{s(t)) is also absolutely continuous. A simple 
computation shows that the energy function Kq ~ U of a is the constant —h for 
almost every t £ [a,/?]- Hence 

V2AjM (7 |[a,b] ) = V2^,7A/(cr) = Aed(cr) " - a). 

But 7 \ [a.b] is a minimizer of Ajm^- So this equation says that a minimizes Ared 
among all absolutely continuous paths in Qh joining (T(a) to cr(/?) in time /3 — a! Let 
^ : [a, /3] — > be a continuous zero angular momentum lift of a to the nonreduced 
space C^. (There are a circle's worth of such paths.) The path ^ is a local minimizer 
of the nonreduced L agrangian action J Ldt among absolutely continuous paths 
joining ^(a) to in time /3 — a. By Marchal's Theorem (see [10] and [4]), ^(t) 
is collision-free for t £ {a, j3) and is a classical solution of Newton's equations. The 
quotient path cr(t) is thus a collision- free solution of the reduced equations with 
energy —h. QED 

To prove theorem [3] we will need certain monotonicity properties for the shape 
potential V. Recall that the shape sphere CP^ endowed with its Fubini-Study 
metric is isometric to (S*^, ds^ /A), where ds^ is the standard metric on the unit 
sphere The isometry is 
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where 

(33) W2 + tW3 = 2^/^1^266 

=/xi|eiP + M2|6p. 

CoUinear configurations correspond to the great circle ^3 = which we call the 
collinear equator. The Northern hemisphere {w^ > 0) corresponds to positively 
oriented triangles. The three binary collision points 6i2,f'23,&3i lie on the equator 
and split it into three arcs denoted Ci, (72,6*3 with the endpoints of Cj being bij 
and bjk- Thus Cj consists of collinear shapes in which qj lies between qi and qk- 

Introduce standard spherical polar coordinates (0, 9) on the shape sphere so 
that one of the binary collision points 6, say b = 612 is the origin {(f) — 0) and so 
that the Northern hemisphere is defined by < 6* < tt. The distance of a point 
from b relative to the standard metric is 4> and is a function of ri2 alone |16j : 
(j) = (/>(ri2);ri2 = ri2{4>) on the shape sphere. Thus setting (f) = const, is the same 
as setting ri2 — const., and defines a geometric circle a = on the shape sphere. 
The collinear equator is the union of the curves 6 — and 9 ~ tt. The arcs Ci 
and C2 are adjacent to b and we can choose coordinates so that Ci is contained in 
the half equator 9 = while C2 contained in the other half equator 9 = n. In the 
case of equal masses, 9 = tt/2 defines the isosceles circle through the collision point 
b = bi2, lying in the Northern hemisphere. 

Lemma 2. Let a be any of the half- circles (p = const centered at 612 and lying in 
the Northern Hemisphere of the shape sphere. Then a can be parameterized so that 
V\a, the restriction of V to a, is strictly convex. Consequently V\a has a unique 
minimum. If that minimum is an endpoint of a then it lies on the 'distant' collinear 
arc C3. Otherwise, the minimizer is an interior minimum and coincides with the 
intersection of cr with the isosceles arc ri3 = r23. 

Remark 1: In all cases, an endpoint of a which lies on arc Ci or C2 is a strict 
local maximum of the restricted V. 

Remark 2: The minimum of l^jo- is an endpoint of a if and only if a does not 
intersect the isosceles arc at an interior point. 

Remark 3: When mi = TO2 the isosceles circle in the Northern hemispher 
(defined by 9 = 7r/2) bisects all the circles a and so the minimum of yjo- always 
occurs at this bisection point and thus is always interior. 

Proof. Use squared length coordinates Si = r"^^ introduced by Lagrange and re- 
introduced to us by Albouy [T]. The Si are subject to the constraints 

> 

and 

(34) 2siS2 + 2S2S3 + 2S3S1 - Si - §2 - S3 > 

which together define a convex cone in the 3-space with coordinates Si. The points 
of this cone faithfully represent congruence classes of triangles, where two trian- 
gles related by reflection are now considered equivalent. The origin of the sec- 
ond inequality is Heron's formula for the signed area A of a triangle: 16A^ = 
2siS2 -I- 2S2S3 -I- 2S3S1 — s\ — S2 — s\. 
In squared length coordinates 

(35) = A^iSi + M2S2 + M3S3 
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(36) U 



'Sl \/S2 -JSS 



A half-circle ct as in the lemma is defined by the two linear constraints — 1, — 
r^2 = c in the Si . Consequently, such a half-circle is represented by a closed interval 



whose endpoints correspond to the coUinear triangles at which the inequality (34) 
becomes an equality. [/ is a strictly convex function of (si, S2, S3) in the region Si > 
and V on the half-circle becomes U (si , S2^ S3) restricted to this closed interval. Now 
a strictly convex function restricted to a convex set remains strictly convex, so 
relative to any affine coordinate on the interval which represents a we see that V 
is strictly convex, as claimed. 

We can parameterize the half-circle cr by 

si(u) = ^(1-Ai2u-Ai3c) 

S2(u) =U 
S3{u) = c 

where u € with u = corresponding to = 0, and u = u+ corresponding 

to 9 = IT. A simple computation gives 

(37) ^V{s,{u),s,{u),s,{u)) = (si(«)-^/2 - s,{u)-'/') , 

which proves that V\cr has an interior critical point if and only if cr cross the isosceles 
curve Si = S2. Convexity implies this crossing point is the unique minimum. 

If the starting point u = u_ of cr is on Ci, then at this point we have 1 between 
2 and 3 so that ri2 +ri3 — r23. Squaring, we see that S2 > si from which it follows 



that the derivative ( 37 ) is negative at u = u_ , and showing that this point is a local 
maximum for ■ A similar argument shows that if an endpoint of a lies on C3 then 
that endpoint is also a local maximum of V^. (If the endpoint of cr, corresponding 
to M = u+ lies on C3 then the derivative ( [37| is positive.) QED 

We prove now continuity properties of the JA/-distance djMiQo, li) introduced 
just before the proof of Lcmmajlj If qo and qi are coUinear we denote by dj]^.j{qo, qi) 
the infimum of the JA/-action among coUinear rectifiable curves in Qh joining qo to 
qi. In a similar way, djj^j{qQ,dQh) denotes the minimum of the JM-action among 
coUinear rectifiable curves in Qii joining qg to the Hill boundary. 

Lemma 3. djM : Qh x Qt ^ -K and djj^j : Ch xCh ^ IR are continuous functions. 

Proof. By the triangle inequality it suffices to prove that for every qo (z Qh or 
Co e Ch the functions q djM{qo,q) 

01 q — ^ d'y jf^j {cq , q)) are continuous in qQ or 
in cq. We give the proof only for q — >■ djM{qo, q)- The proof for q — > d'^jj^j(qo, q) is 
similar. 

If go is a non-collision shape, not lying on the Hill boundary the djMiqoi <?) is a 
regular Riemannian distance in a neighborhood of qo, so the result is classical. If 
qo lies on the Hill boundary, the function is a pseudo-distance: d{qo,q) = with 
Qa 1 can happen), but the classical result goes through with no changes. So 
assume now that qo is a collision point. We consider two case. 
1st case : qo is the total collision. Let us choose any non-collision coUinear 
configuration on the shape sphere s ^ S^. If q — {r,s) £ Qh, and 7 is the path 
obtained by joining the total collision to g = (r, s) along the ray through q, and 
then joining g to g along a segment of geodesic of {r} x S^, a simple computation 
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shows the existence of a constant C > 0, independent of g, such that Ajm{'^) < 
Cy/r + 0{r^/r), hence djMiqo-.q) ^ as q ^ go- 

2nd case : go a partial collision. Assume that go is on the coUision ray ri2 = 0. 
Introduce spherical coordinates (0, 9) as before centered at the binary cohision 612, 
so that (f> — corresponds to 612, and 9 — or tt corresponds to the coUinear circle. 
If ro is the radial coordinate of go and gi G Qh is a point with radial coordinate ri 
and spherical coordinates {(j>i, 9i), let us choose a g (0, tt) and define the following 
three paths in spherical coordinates 



71 
72 

73 



r(<)=ro, m^t^, 0{t)^9i, i G [0, 1] 

r{t) = {l-t)ro + tri, (b{t) = 4>, 9{t) ^ 9^, t G [0, 1] 

r(<)=ri, = + 9{t) = 9^, tG[0,l]. 



Let 7 be the path obtained by pasting together 71, 72 and 73. This 7 joins go to 
gi . By ([6]) and ( 33 1 we get 



V{<j^,9)^l+o(h, 

where c > depends only on the masses. If rf)/2 < ri < 3ro/2, a simple computa- 
tion shows that 

Ajm{ii) + 

hence chosing = |ri — ro| we get 

djM{qo,qi) < AjMh) < C\/\rx - rol + o(Vki - ?'o|), 
where C > depends only on ro, hence djM{qo, gi) — > as gi — > go. QED 

Proof of Theorem |3} Let go G C/i be a coUinear configuration. By Proposition 
15 there exists a JAf -minimizer 7 : [0,r] — > among paths starting in go and 
ending on the Hill boundary dQh- The JM- action of a path on the Hill boundary 
is zero, therefore if ^o is the first time where 7(to) G dQh, then ^{t) G dQh for 
all t G [to,T]. We may cut out the arc 7|[to,T] from 7, or equivalently, assume 
that 7(i) ^ dQh for all t < T. For every t G (0, T) the restriction 7 | [o^tj is a 
fixed end point JAf-minimizer. By Lemma [l] 7 is a classical brake solution, and is 
collision- free for t G (0,T]. 

Let us prove now that if 7 is not coUinear at all times, it has a unique syzygy : 
the point go. Indeed, assume for the sake of contradiction that 7(t) is coUinear for 
some < r < T. Recall that the set of initial conditions tangent to the coUinear 
submanifold yield coUinear curves. CoUinear brake points (with zero velocity) are 
such initial conditions. Consequently, r = T is impossible, for if 7(T) G dQh 
were coUinear, all of 7 would be coUinear. So t < T. By the same reasoning, 7 
cannot be tangent to the coUinear manifold at t < T. Consider the path obtained 
by keeping 7 | [o^t-j as it is and reflecting 7 | j-j with respect to the syzygy plane. 
This new curve has the same JM-action as 7, and joins go to dQh, so it is a JM- 
minimizer. But this new path is not differentiable at t = r, contradicting the fact 
that minimizers are solutions, hence smooth. We can now assume without loss of 
generality that 7 is all the time in the half-space of positive oriented shape. 
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If qo is the total collision wc will now show that 7 is the Lagrange homothety 
solution. Write 7(t) = {r{t),s{t)) G Qh <Z M+ x S^. Denote by sl e the 
Lagrange shape in the Northern hemisphere. We want to show s{t) = sl for all t. It 
is well known that the Lagrange shape is the unique minimum of the shape potential 
V. Let be the number such that :pj^V{sl) = h. Suppose that s{T) 7^ s^. Then 
r{T) > TL since V{s{T)) = r{T)h > V{sl) = rLh. Since r(0) = there will be a 
smallest number t = T.^, in the interval [0,r] such that r{t) ~ rj^. Then r{t) < rj^ 
for t < T^. The path ^L{t) = {I'i't), sl), < ^ < lies in Qh, joins go to dQh and 
satisfies C/(7L(i)) ^ U{'^{t) on this interval (with equality if and only if s{t) = sl)- 
The kinetic energy of this new path is pointwise the same or less than that of the 
old path, over their common domain. Consequently ^ja/(7l) < Ajm{^), with 
equality if and only if s{t) — for all t G [0,T] (and = T) as desired. 

Assume now go is a double collision, say on the collision ray ri2 = 0. We prove 
that our minimizer 7 is not coUinear. Assume, for the sake of contradiction that 
^{t) is collinear at all time t G [0,r] . Since 7 is collision free except for q^, 7 
is contained in an angular sector of the syzygy plane between the collision ray 
ri2 = and one of the two other collision rays. Say it lies in the sector Ci between 
— and — 0. Introduce spherical coordinates (</>, 9) on the shape sphere 
as in Lemma[2] Then 7(t) has coordinates {r{t),(j){t),9{t)) — {r{t), (/){t),0), where 
we have arranged the coordinates so that the sector within which 7 lies, sector Ci, 
is given hy 9 = 0. To re-iterate, at every moment 7's spherical projection lies in 
the arc Ci on the collinear equator. By Lemma [2] there exists S G (0, tt) such that 
V{(l){t),6) < V{(l){t),0) for all t G [0,T]. Let ts be the first time in our interval 
such that V{4>{t),S)/r{t) = h. (There is such a time by an argument similar to the 
one above when qo was triple collision.) Define the new path ^sit) = {r{t), 4>{t), (5), 
t G [0; Ts\. Geometrically, 75 is obtained by rotating ^{t) by an angle 5 around the 
collision ray ri2 = and keeping only that part of it which stays inside Qh- Such 
a rotation is an isometry for the metric ([T]), but strictly decreases U. Therefore 
^jm{is) < -^juil)- But 7 is a minimizer, so we have a contradiction. 

Using the notations of Lemma [3] we have proved that if go is a collision, we have 
djM{qQ,dQh) < d'}j^,j{qo,dQh)- 

Assume now that go G Ch is a collision configuration. We show the existence of 
a neighborhood O (go ) of go in C'h such that for every go G O (go ) , no JM-minimizer 
from go to dQh is collinear. Indeed, let us assume, for the sake of contradiction, 
the existence of a sequence (gn)neN in Ch, converging to go, with corresponding 
collinear minimizers. Then djM{qn, dQh) — djAiiln, dQh) for all n. By the triangle 
inequality we have 

\djM{qn,dQh) - djM{qo,dQh)\ < djM{qn,qa), 
\djMiln,dQh) - d';j^.j{q^,dQh)\ < d'}^j{qn,qa)- 

Applying Lemma [s] we get rfjM (go, dQh) = dj]^j{qo, dQh), which is a contradiction. 
The neighborhood U of of the collision locus described in the Theorem can be taken 
to be the union of the C(go) as g^" varies over the collision locus. This finishes the 
proof of Theorem [3] 

Proof of Theorem [4| In order to prove part (a) of the Theorem, assume mi — 
m2 and let go be a configuration on the collision ray ri2 — 0. Let 7 : [0, T] Qh 
be a JM- minimizer from go to dQh- If go is the triple collision point, by Theorem 
[3] we know that 7 is the Lagrange brake solution, which is an isosceles brake orbit. 
If go is a double collision point, then express 7 in terms of the spherical coordinates 
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system {r,(j>,9) on the shape space S^, centered at the double coUision 612 as m 
Lemma [2] and in the proof of Theorem [sj We can assume that j{t) hes in the 
Northern hemisphere on the upper half-space < 9 < n. Because mi = 7712, the 
isosceles curve ri3 = r23 in the Northern hemisphere corresonds to the great circle 
{6 ~ 7r/2}. By Lemma [2j for fixed 0, the absolute minimum of 9 1-^ V{(j),9) is 
achieved at the isosceles curve 9 = tt/2. Let < t < T be the first time such 
that V{(l){t),'K/2)/r{t) = h. Then the path j{t) = (r(t), 7r/2), t € [0,t], 
joins go to dQh, has U{'j{t)) < U{'^{t)) with equality if and only if 7 = 7 on 
its domain, and has kinetic energy less than or equal to that of 7 on its domain. 
Thus Ajm{'j) < Ajm{i), with equality if and only if ^{t) is isosceles at all times 
t G [0, T]. Since 7 is a minimizer, it is necessarily isosceles. This proves part (a). 

To prove part (b), assume < r23 at the starting point q^. The isosceles set 
fiz — ^23 is invariant by the flow. If we assume, for the sake of contradiction, that 
7 : [0, T] — >■ Qh intersects ri3 — r23 at some time < r < T, then it is necessarily 
transverse to ri3 = r23, otherwise ^{t) would be isosceles at all times. The reflection 
with respect to the isosceles plane ri3 = r23 is a symmetry for the kinetc energy 
and for the potentiel U , therefore the path 7 obtained by keeping 7 | [q.t] as it is, 
and reflecting 7 |[r,T] about the isosceles plane ri3 = r23, is still a JAf -minimizer. 
But 7 is not differentiable ai t = t. This is a contradiction. 

The proof of part (c) is obtained applying (b) to the three isosceles plane = 
r^k■ QED. 



Proof of Proposition |14[ Seifert's Tubular Neighborhood Theorem. 

Our construction is the same as Seifert's 22J. The one real difference is that our 
Hill boundary dQh is not compact, whereas his is compact. 

We will construct Seifert's neighborhood in the non-reduced configuration space 



and then project it via the quotient map tt of eq ( 10 ) to Qh to arrive at the 
reduced Seifert neighborhood. We abuse notation by using the same symbol Qh 
for the reduced and non-reduced Hill regions and similarly using the same symbol 
Sh C Qh for the reduced and non-reduced Seifert neighborhoods to be constructed. 

The configuration space is endowed with the mass inner product < , > which 
is the real part of the Hermitian inner product of cq. ([t]). The equations there are 

(38) 9-VC/(g), 

where the gradient of U is calculated with respect to the mass inner product. The 
Jacobi-Maupertuis action of a curve 7 : [^0,^1] — ^ Qh is defined by 



Ajm{i) = / VKV'^{U-h)dt, 



to 

where K = ||7|P/2. 

As a first step we construct a map F : dQh x [0, S] — )■ Qh which is an analytic 
diffeomorphism onto its image and is such that the curves t 1— F{x,t'^) are the 
brake solutions starting from x a.t t = 0. Given a < h < f3, define the open set 

-Da,/3 = {a; e C^, a< U{x) < P}. 

If a; e Da.p, then the smallest of the mutual distances rij of the triangle defined 
by X is bounded below by a constant depending only on (3 and the masses. We 



denote by Qxit) the solution to (38) with initial conditions qx{0) — x and qxiO) — 0. 



(Note: these solutions need not have energy h.) Choose positive numbers a and b 
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such that a<a<h<l3<b. By classical results on differential equations, there 
exists T > such that for every x £ D^^p, the solution qx{t) is well defined for 
t G [— T, T] and satisfies Qxit) € Da.t- The map (x,t) i— qx{t) is analytic, even in 
t (i.e. qx{—t) = qx{t)), and its derivatives up to the second order are uniformly 



bounded on -Dq,^ x [— T, T]. By equations of motion (38) 

qx{t)=x+^^t'+0{t^). 

Let us set Ti = \/T. Since qx{t) is even in t, the map F{x,t) = qxi^^/r), {x,t) G 
Da,i3 X [0, Ti] is still analytic, so it can be extended to negative values of r, giving 
an analytic map F : D^./^ x [^Ti,Ti] — )■ Da^b satisfying 

(39) F{x,T)^X+'^^T + f{x,T), 

where the f{x, r) — 0{t^), uniformly for x £ 
Define 

G : X [-Ti,ri] ^ Da^b x K, G(x, r) = (F(a;, r), C/(x)). 
By (39), the differential of G at a point (x,t = 0) gives 

-DG(^^o)('5a;,(5r) = (fe H ^St,< VU{x),Sx >). 

An easy computation shows that DG(xfl) is invertible and || (_DG(2,.o))^^ II is bounded 
as a: S D^^p. Moreover, since all derivatives of G up to the second order are uni- 
formly bounded on D^ j^ x [— Ti,Ti], we can apply a strong version of the inverse 
function theorem and find a (5 > such that for all x e -Dq,^, the map G defines 
a diffeomorphism from B{x,5) x [—(5,(5] into its image, where B{x,d) denotes the 
closed ball centered in x with radius 6. If we take the restriction to U{x) = h and 
define BgQ^{x,S) = B{x,S) r\ dQh we find that F defines an analytic diffeomor- 
phism from Bqq^ {x, S) x [—S, 5] into its image. Let us prove now that by decreasing 
sufficiently 6, the restriction of F to dQh x [—S, 5] is an analytic diffeomorphism onto 
its image. Indeed, we have proven that at every point it is a local diffeomorphism. 
Assume, for the sake of contradiction, there exist two sequence (a;„,r„)„gN and 
(x^,r4)„gN satisfying (a;„,T„) ^ (a:J^,T,'j) such that t„ — >• 0, — >■ as n — ^ +00 



and F{xn,Tn) = F{x'„,T^) for all n e N. By (39) and uniform boundedness of VC/ 
on dQh, we have ||x„ — a;^|| — )■ as n +00, therefore x'„ S BgQ^{xn,S) if n is 
sufficiently great. This contradicts the fact that -F is a diffeomorphism on every 
BaQ^{x,S) x [-6,6]. 

The second step consists in defining a diffeomorphism <& on dQh x [0, 6] of the form 
$(x,?/) = F{x, /3{x,y)) where the scalar function /3 allows us to "straighten out" 
the JM action. Let S{x,t) denotes the JAf-length of the extremal F{x, s)s£[o,t]- 



From ( 39 ) we compute 



(40) S{x,r)=T'^'(^^\\VU{x)\f+gix,T) 

where g{x,T) is analytic, and g{x,T) — 0{t), uniformly on a; G dQh- The map 
H{x,t) = (x, 5(0;, r)^/^) is analytic in dQh x [0,(5], and we can find e > such 
that H maps a neighborhood of dQh x {0} (in dQh x [0, 6]) diffeomorphically into 
dQh X [0,e]. Then H'^ has the form {x,y) ^ ix,l3{x,y)). We set $ = o H-^ 
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where the domain of H ^ is dQh x [0, e]. It is clear now that for every x E dQh, 
and for every S e (0, e), the curve {^{x,y))y^[Q S] is a segment of brake orbit and 
its JM-action is S^/^. 

The Jacobi metric in Qh C is given by 

g,,ix){v,v)^2{U{x)-h){v,v) 

We now show that brake solutions {^{x, y))ye[o,e] are g/i-orthogonal to hypersurfaces 
^{dQh X {y}). The argument is the standard one used to prove the Gauss lemma 
in Riemannian geometry. We are to show that gh{dy^(x, y), dx^{x, y)£,) ~ for all ^ 
tangent to the hypersurface. For this purpose, set f{y; x) = gh{dy^{x, y), dx^{x, y)£,) 
for fixed Because the lengths of the curves y — )■ ^{x, y) up to the value y^ are all 
y^^ we have that ghidy^ix,y),dy^{x,y)) — |y, independent of x. Differentiating 
this identity in the ^-direction tangent to the hypersurface and commuting deriva- 
tives we get gh{dy^{x,y),\/ydx^{x,y)£,) = where V denotes the g/j-Levi-Civita 
connection. Now, because the y-curves are reparameterized geodesies with lengths 
only depending on y we have that \7ydy^{x,y) = fi{y)dy^{x,y) for some function 
fi{y). Differentiate with respect to y, and use metric compatibility and the 

previous equations to derive the linear differential equation dyf = ^J'{y)f{y',x) with 
a; as a parameter. But /(O; a;) = so the initial condition of this differential equa- 
tion is zero, from which it follows that / is identically 0. (Indeed /i is the second 
derivative of y^^^ with respect to y, which is singular at y = 0, but the analysis 
goes through.) 

An arbitrary tangent vector v £ TqQh a,t q — $(2;, y) can be written 

V = dx^{x, y)^ + dy<^{x, y)A, 

where (^, A) G T^dQh x M. By the previous orthogonahty discussion and the fact 
that y^/^ is the arclength of ($(a;, y))j,g[o,e] have for such v: 

(41) gk{v,v) = {U- h)Kn(x,y){i,S) + ^A^. 

where Kfi{x,y) is a positive definite quadratic form on TxdQh- If 7 : [^Oj^i] ~^ Qh 
is any rectifiable curve joining a point q — ^{x, d) (with < (5 < e) to dQh, and if / 
is the set of times t E [io,ii] such that 7(t) E Sf^ = ^{dQh, [0, e]), by (41) we have 

(42) AjMh) - V9k{j{t),i{t))dt > ^ ^^|y(i)Mt > S'/\ 

where we term {x{t),y{t)) = ^~^{j{t)) for t E I. Moreover, we have equality in 
(42 1 if and only if, up to an arc contained in dQh, the curve 7 is a reparametrization 
of {^{x,y))y^[Q^s]- The Euchdean length of (^(a^, y))yG[o,<5] is given by the integral 
/o \\dy^{x,y)dy\\ dy. By (39), (40) and by definition of $ we have 

Since \7U{x) is uniformly bounded on dQh, there exists a strictly positive constant 
M, independent of x and d as long as 5 < e, such that 

\\dy^x,y)dy\\dy<MS. 
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Let us show now there exists a > such that for every 5 G (0, e) we have U > h+ad 
on Qh\Sl. By we have 



(44) U^x, y)) = h + " ^ y^!, y + 0{y'), 

for {x,y) € dQh x [0, e]. Since ||V[/(a;)|| is bounded and uniformly bounded below 
below by a positive constant, we can find two constant < a < /3 such that 

h + ay < y)) <h + Py 

for every {x,y) G dQh x [0, e]. Assume now, for the sake of contradiction, the 
exisence of a point q e Qh \ such that U{q) < h + a6. Since level set of U are 
connected by arc, there would exist a point q' on ^{dQh, S) such that U{q') < h+a6, 



and this is in contradiction with ( 44 1 



5. Periodic Brake Orbits 

The goal of this section is to prove Theorem [2] about the existence of simple, 
periodic brake orbits for the isosceles three-body problem. Other examples of pe- 
riodic, isosceles brake orbits, more complicated than the one described here, are 
given in pS] , 

Let p be a brake initial condition other than the Lagrange homothetic one, po- 
According the results of section |3.H p can be followed forward to meet the syzygy 
submanifold, C'h, i.e., the submanifold of the energy manifold with collinear shapes. 
There is a natural reflection symmetry of the energy manifold through Ch, obtained 
by reflecting the shape variables {x,y) and their velocities {x',y') while leaving the 
size variables (r, unchanged. Call this reflection map R. Then i? is a symmetry 
of the differential equation if one also reverses time. If the orbit, 7(t), of p meets 
Ch orthogonally after time Ti , then reflecting the orbit segment and reversing time 
gives the continuation of 7 to the time interval [Ti,2Ti] and 7(2r) = R{p), the 
reflection of p. Now it follows from symmetry during the interval [2Ti,4Ti] the 
orbit retraces its path, returning to p. So p determines a periodic brake orbit of 
period AT. 

If j{t) does not meet Ch orthogonally and if jiT) is not in the local stable 
manifold of the Lagrange triple collision, then it can be followed to a second syzygy, 
say at time T2 > Ti. If this crossing is orthogonal one obtains a periodic brake 
orbit of period by reflection. 

We will see that such a second-syzygy periodic brake orbit exists in the isosce- 
les subsystem of the three-body problem, at least for certain choices of the mass 
parameters. It is not known whether a first-syzygy brake orbit exists. Numerical 
experiments suggest that no such orbits exist in the equal mass case. 

5.1. The Isosceles Three-Body Problem. Assume that two masses are equal, 
say mi = m,2 = 1. Then there is an invariant submanifold of the three body problem 
such that the shape remains an isosceles triangle with ruj, on the symmetry axis for 
all time. Up to rotation, this isosceles subsystem is obtained by making the Jacobi 
variables ^1,^1 real and £,2,^,2 imaginary. 

As in section [2] separate size and shape variables will be used. In addition, double 
collisions will have to be explicitly regularized. A convenient way to do this is to 
define an angular variable which gives a multiple cover of the isosceles shape space 
and which is locally a branched double cover near the binary collisions, as in the 
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familiar Levi-Civita rcgularization. Using the projective Jacobi variables [^i,i^2], 
this parametrization of the isosceles shapes can be accomplished by setting 

6-^cos2(0) ^2-^2sin(0). 



The isosceles binary collision shape corresponds to = or 6* = ±-| mod 2tt. 
Note that with the assumptions about the masses, one has 



1 2m 



3 



2 2 + mg 

Substitution gives a reduced Lagrangian 

Lredir,r,0,0)^Ko + -V{9) 
r 



where 
and 



2r2 cos^ e ^2 



(1 + sin' 



2 , 



^ cos2 9 + sin^ 6')2 + ^ sin^ ( 



3 

Introducing a change of time scale ' = r 2 cos ' leads to the following system of 
differential equations: 



(45) 



it;2 cos2 e + lw^{l + sin2 Of - W^e) 



\w{l + sm" 



l2 0)2 

= W'{e) ~ lvwcos^9 + smecose{2r + v^ - ^w^il + sm^ 6)) 

where W{6) — cos^ dV{9). The energy conservation equation is 

(46) iw2 cos2 9 + |u;2(l + sin2 9)^ ~ W{9) = cos2 9 

where the energy has been fixed at —1. 

Since W{0) is an analytic function, the binary collisions have been regular- 
ized. Moreover the triple collision singularity has been blown up into an invari- 
ant manifold at {r = 0} as before. The isosceles Hill's region for energy —1 is 
Qi — {{r,9) : < r < V{9)}. This is shown in figure [6] for the equal mass case 
mg = 1. 

Using these coordinates, syzygies occur at the Euler shape {9 = mod 27r) and at 
binary collision (6* = ±^ mod 27r). These are indicated by the bold vertical lines in 
figure |6] The reflect ion through syzygy amounts to reflecting the position variables 
{r,9) through these vertical lines while taking the velocity (w,w) to {v,—w) and 
reversing time. To get a symmetric periodic orbit, one needs to reach syzygy with 
V = 0. Since r' = vr cos^ 9 this is equivalent to orthogonality at the Euler vertical 
lines 9 = mod 2tt, but every orbit crosses the lines 9 = zL^ mod 27r orthogonally. 
A numerically computed periodic brake orbit is shown in the figure. It begins on 
the zero velocity curve, crosses the syzygy line at 9 = with v < (not apparent 
from the figure), then continues to its second syzygy at 9 = where it crosses 
orthogonally with v = 0. The rest of the orbit is obtained by symmetry. The goal 
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Figure 6. Isosceles Hill's region for the equal mass three-body 
problem using coordinates {0,r). The zero velocity curve is the 
top boundary curve. The syzygy configurations are represented by 
the thick vertical lines (61 = is the coUinear central configuration 
and 9 = ±^ are binary collision shapes) . A numerically computed 
periodic brake orbit is also shown. 



of this section is to prove the existence of such an orbit for certain choices of the 
parameter m^. 

The idea of the proof can be described briefly as follows. Consider a curve, Z, 
of brake initial conditions whose shapes vary from equilateral to binary collision. 
In figure |6j the equilateral shapes (minima of the shape potential) are shown with 
dashed vertical lines. Let 6* be the equiliateral shape in the interval [0, Other 
equilateral shape occur at the points ±9* + kir where k is an integer. The curve Z 
will consist of the brake initial conditions with 9 £ [9* — tt,— (in figure [61 this 
is the part of the top boundary curve between the left-most dashed and left^most 
bold vertical lines). It will be shown that as the initial condition p E Z varies, 
there is at least one point which can be followed to meet the vertical line 9 — 
orthogonally. 

The proof will use a geometrical argument in the three-dimensional energy man- 
ifold, Pi = {(r, 9,v,w) : r > 0, H — —1}. This manifold can be visualized through 
its projection to (r, u, 0)-space, which is given by the inequality r + ^v^ < V{9). 
This projection is shown in figure [7) On the top surface of the projection, w — 0. 
The full energy manifold can be viewed as two copies of this projection (one with 
■w >0 and one with w < 0) glued together along this top surface. The desired orbit 
has the property that the shape angle 9{t) will increase monotonically from 9{0) to 
9{T2) — where T2 denotes the second-syzygy time described above. Therefore it 
suffices to consider the part of the energy manifold with 9' — jw{l + sin^ 9)^ > 0. 
On this half of the energy manifold, one can solve (46) uniquely for w{r,v,9). 
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Figure 7. Projection of the w > part of the energy manifold 
in coordinates 9 (width), v (depth), r (height), The top surface is 
{w ~ 0}; the floor is part of the coflision manifold {r — 0}. A 
numerically computed periodic brake orbit is also shown, passing 
through the walls of the regions Ri , Ru , Rm as described in the 
proof. 



The relevant part of the energy manifold will be divided into three regions: 
Ri=Pin{9e [61* -7r,-f],w > 0} 

Rii = Pin{e& [-f ,-r],u) > 0} 
Riii = Pin{9 e [-r,o],w>o}. 

The planes 9 — 6* — n and 9 = will be called the left and right walls of Rj 
respectively with similar definitions for the other regions (the walls are the four 
vertical planes in figure [t]). 

These regions contain certain restpoints on the collision manifold which will now 
be described. Let L± be the Lagrange restpoints in Rj at (r, 9, v, w) = (0, 9* — 
TT, ±i;*,0) where v* = y/2V(9*). These are connected by the Lagrange homothctic 
orbit which is the curve of intersection of the left wall of Ri with the boundary 
surface {w — 0}. Similarly, there are Lagrange restpoints L'_j. at {r,9,v,w) — 
(0, —9*,ztv*, 0) and a corresponding homothetic orbit in the left wall of Rm which 
is also the right wall of Ru. Finally there are Eulerian restpoints E± at (r, 9, v, w) = 
(0, 0, ±^2F(0), 0) and an Eulerian homothetic orbit in the right wall of Rm- 
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Let Z be the part of the zero velocity curve in Rj. The proof will follow a subset 
of Z forward under the flow through these three regions to obtain a curve in the 
right wall of i?///, i.e., in the syzygy set at 6* = 0. It will be shown that this final 
curve crosses the plane u = and the point of crossing will determine the required 
periodic brake orbit. The next lemma shows how the flow can be used to carry 
orbits across the various regions. 

Lemma 4. Regions Ri and Rm are positively invariant sets for the flow while 
Rii is negatively invariant. With the exception of the Lagrange homothetic orbits, 
orbits cross these regions as follows: any orbit beginning in the left wall of region 
Ri, crosses the region and exits at the right wall. The same hold for Rm except 
for orbits in the stable manifold of ( which is contained in the triple collision 
manifold {r — 0}). 

Similarly, except for the Lagrange orbit, any backward-time orbit beginning in 
the right wall of Rn can be followed back to the left wall. Finally, forward orbits 
beginning in the left wall of Rjj either leave Rn through the right wall, leave Rn 
through the top surface {w ~ 0} or converge to one of the Lagrange restpoints L'^ 
in the right wall as s ~> od. 

Proof. By definition, an orbit in any of the three regions satisfies 9' > so 0{s) 
is non-decreasing. Referring to figure [7j the lower boundary surface {r — 0} is 
invariant. On the upper boundary surface, w = and 

w' = w'{e) + smecose{2r + v^) ^ w'{t) + 2smecosewie) = cos^ev'{e). 

Now Ri,Riii are are defined by ^-intervals where V'{9) > so for these regions, 
w' > on the top boundary. This proves positive invariance and also that it is only 
possible to leave through the right wall. Similarly, in Rjj we have w' < on the 
top wall so the region is negatively invariant and backward orbits can only leave 
through the left wall. 

Assume for the sake of contradiction, that an orbit other than the Lagrange ho- 
mothetic orbit remains in Ri for all time s > 0. Then 9(s) converges monotonically 
to some limit 9^0 G {9* — vr, — ^] (it cannot be ^* — tt since all points in this plane 
but not on the Lagrange orbit have w > so they initially move to the right) . The 
omega limit set is either empty or else it must be a nonempty invariant subset of 
{9 — 9oo}- However, these planes don't contain any nontrivial invariant sets. So 
the omega limit set must be empty which is only possible if 0oo = ^-i^d if the 
orbit leaves every compact subset of the energy manifold. To show that this is 
impossible, let A = ^/2r + v'^ so that the energy equation becomes 

cos^ 9 + lw^{l + sin^ 9)^ = W{9). 



We will show that A(s) remains bounded to complete the argument. From (45 1 we 
find AA' = ^vw {l+sm'^9f. Since 9{s) is increasing, we can reparametrize by 9 
to get ^ = 1^ and since < A we have 

m i§i<i- 

As w{s) is bounded by the energy relation, we get a bound for |^|. It follows that 
A(s) is bounded along the part of the orbit in Rj as required. 

The same argument applies to backward-time orbits in Ru, showing that they 
must reach the right wall. The argument for forward orbits in Rm is easier since 
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this region is compact. The omega hmit set of an orbit remaining in Rm for all 
time would have to be a nonempty invariant set in a plane — € (~fiO]- 
The only invariant sets are the Eulerian restpoints and homothetic orbit and so 
the omega limit set would have to be one of the restpoints. However, the stable 
manifold of is just the restpoint itself and the homothetic orbit, so the orbit 
must be in the stable manifold of which, it so happens, is contained in the 
collision manifold. QED 

In addition to this lemma about region-crossing, it will be necessary to use some 



facts about the flow on the isosceles triple collision manifold. Setting r = in (45) 



gives the dynamics on the triple collision manifold. The energy equation ( 46 1 gives 



cos^ 



' + \w^{l + sin^ Of - W{6) = 0. 



Using this to eliminate w gives a flow on part of the (0, v) plane satisfying < 
2V{9). This is shown in figure [s] for the case ma = 1. An important property of 
the flow is that it is gradient-like with respect to the variable v. Indeed, using the 
energy equation gives v' = ^^^(l -I- s'm^ 9f > and it can be shown that v{s) 
is strictly increasing except at the restpoints. The restpoints L±,L'j. are saddle 
points. 

Certain properties of their stable and unstable manifolds will be used in the 
proof. Let 7,7' denote the branches of W^{L-), VF"(L'_) in {w > 0} (bold lines in 
figure [8| . The key properties needed to complete the existence proof refer to the 
intersections of these branches with the syzygy lines at 6* = — ^ and 6 = 0. We 
require that 7 remains in {w > 0} at least until it crosses these two syzygy lines and 
that the intersection points are {6,v) = (— where t;i < and {0,v) = (0, U2) 
with V2 > 0. Furthermore we require that 7' remains in {w > 0} at least until it 
crosses {9 = 0} at a point (0, w) = (0,W3) with V3 < 0. Call a mass parameter TO3 
admissible if these hypotheses hold. The figure indicates that 7713 = 1 is admissible 
and the next lemma guarantees that this is so. The proof will be given later. 

Lemma 5. There is a nonempty open set of admissible masses which including the 
equal mass value = 1. 

These lemmas can be applied to follow a subset of Z through Rj to a curve 
Zj in the right wall. The left endpoint of Z is the brake initial condition with 
equilateral shape, 9 = n — 0* . Call this point pq. The corresponding orbit is the 
Lagrange homothetic orbit, so it remains in Rj for all i > and converges to the 
restpoint L_ . It follows from lemma |4] that Z \ po can be followed through i?/ 
to the right wall. It will be important to understand the image curve Zi whose 
three-dimensional projection lies in the half-plane {{r,9,v) : r > 0,9 = — f}. Use 
(w, r) as coordinates in this half-plane and let pj be the point (w, r) = (vi, 0) where 
the branch 7 of W'^{L^) meets the right wall. 

Lemma 6. Assume 771,3 admissible. Then the image curve Zj of Z \ pq is a 
continuous open arc in the right wall of Rj . One end converges to pj = (ui, 0) and 
at the other end |(7', u)| — > 00. 

Proof. The initial curve Z \po is a continuous open arc in the zero velocity curve 
from Po to {r,9,v,w) = (cx), — |, 0, 0). The flow across the right wall of Rj is 
transverse since w > there. So Zj is a continuous open arc. Initial conditions 
near po on Z will follow the Lagrange homothetic orbit to a neighborhood of the 
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Figure 8. Flow on the w > part of the colhsion manifold in the 
equal mass case (the "floor" in figure]?]). The crucial branches 7, 7' 
of H^"(L„), W"(L'„) are shown. 7 intersects the hne 9 = with 
V < and 6 = with v > and 7' intersects 9 = with t; < 0. 
Two shorter branches of stable manifolds which play a role in the 
proof are also shown. 



restpoint L_ and then follow the branch 7 of the unstable manifold to meet the 
right wall near the point pi. 

Next consider an initial point near the other end of t he c urve. Initially the 
quantity A = \j2r + is large and « — ^ . It follows from ( 47 ) that A is still large 

QED 



is large and « — ^. 
when the orbit meets the right wall. 



Now part of the curve Zi will be followed across region i?//. Unfortunately, this 
region is not positively invariant (solutions can leave by w becoming negative, i.e., 
Q begins to decrease). It turns out, however, that part of the curve Zi is trapped 
inside Ru by an invariant surface, namely, the stable manifold W^(L'_). From the 
linearization at the restpoint, it follows that W*(L'_) has dimension 2. One of the 
orbits in the stable manifold is the Lagrange homothetic orbit which lies in the right 
wall of Rii connectin to L'_. It is known that W^(L'j^ and W^iL'S) intersect 
transversely along this orbit [21]. W^iL'J) also contains two orbits in the collision 
manifold, one of which lies in region Ru. The one-dimensional manifold W^iL'^ 
also contains an orbit in Ru . These two branches of stable manifolds are shown as 
dashed curves in figure [S] 

Consider the "quadrant" of the surface W^(L'S) in region Ru- One edge is 
the orbit in the collision manifold just described and the other is the Lagrange 
homothetic orbit in the right wall. It follows from lemma [4] that with the exception 
of the homothetic orbit itself, orbits in this quadrant can be followed backward 
under the flow to reach the left wall of the region. The intersection of the surface 
and the wall will be a curve. One endpoint of the curve arises from the branch 
of W^{L!_) in the collision manifold. Since v is decreasing for backward orbits 
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in the collision manifold, this endpoint will be of the form {v,r) = (vcO) with 
vq < —V* < vi < 0. To find the other endpoint, note that backward orbits in the 
quadrant near the homothetic orbit will follow the homothetic orbit back near the 
restpoint L'^ and then follow the branch of W'^{L'_^) in Rjj. This branch of stable 
manifold is related by symmetry to the branch 7 of W'^{L^). Hence it meets the 
left wall of Rjj at (u,r) = (— ui,0), —vi > 0. 

So the surface VK*(L'_) intersects the left wall of Rjj in a curve connecting the 
two points (tiQ, 0), {—vi, 0) in the collision manifold but otherwise lying in {r > 0}. 
It follows from lemma [6] that Zj crosses this curve. Let Zj denote the part of Zj 
below the stable manifold from {v,r) = {—vi,0) to the first intersection, call it qj, 
with W{L'_). Points in Z'j \ qj can be followed forward across region Rjj to its 
right wall. 

Indeed if p G Zj \ qj, the forward orbit of p remains in the part of Rjj below 
the invariant manifold W'^{L'_). By lemma |4| it either reaches the right wall or 
converges to L'^.. By construction, the only point of Z'j in W^*(L'j_) is qj. Points of 
Zj near qj will reach the right wall of Rjj near L'_ . The other endpoint of Z'j is the 
point pj with {v,r) = (wi,0) in the branch of iy(L_) described in the definition 
of admissible mass. This point continues to follow that branch to its intersection 
point with the right wall. Thus the image of Zj \ qi is an open arc Zjj in the right 
wall of Rji connecting L'_ to 7 and otherwise lying in {r > 0}. 

Since Zjj does not intersect the Lagrange homothetic orbit or the stable manifold 
of lemma [4] shows that it can be followed across region Rjjj to the right wall 
to form an open arc Zjjj. The endpoints of Zjjj are at the point = (^'2:0) 

and (w3, 0) where the branches of unstable manifolds 7 and 7' cross the wall. Since 
ms is admissible, we have ^3 < < W2, so the is at least one point of Zjjj with 
w = 0. Thus there is a point p of the zero velocity curve Z which can be followed 
through all three regions to reach 9 = Q with w = (the curve in figure 

[7|- Using 

the refiection symmetries we get a periodic brake orbit, completing the proof of 
Theorem m 



5.2. Proof of Lemma [5l The behavior of the stable and unstable manifolds on 
the isosceles triple collision manifold as the mass ratio 7713 varies has been studied 
using a combination of analytical and numerical methods by Simo ([^). His results 
imply that the admissible masses, 7773, are those in the open interval (0, 2.6620). In 
this section we will only prove that = 1 is admissible. The admissible masses 
form an open set so they will then include some open interval around 1. Using 
(0, v) as coordinates on the collision manifold and using the energy equation with 



r = in ( 45 1 gives 

„' = 1^2(1 + sin2 ef 
0' ^ liv{l + sm^0f. 
Any orbit segment with w > can be parametrized by 6 and we have 



dv _ ^ _ 2y/2W{9) ^v'^cos'^e __ 2| cos 01^2^(6*) 
" 2^ ^ 1 + sin^ e ~ 1 + sin^ 6 



Let 7773 = 1. We need to follow the w; > branches 7,7' of the unstable man- 



ifolds W"{L-),W"{L'__). These are solutions of (48) beginning at L-,L'_ whose 



coordinates are {e,v) = (tt - 6**, -w*), (-6**, -w*) where v* = y/W{¥) = ^6. 
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The branch 7' is short and easy to un derstand using crude estimates. Let v{9) 
denote the corresponding solution of (48). We have v{—6*) = —v* = — \/6 and we 
need to show that v{0) < 0. In the interval [-9*,0] we have V{e) < V{0) = 

From (Esl we get ^ < y/2V{0) and 

v{0) < -~v* + ^/2V{0)e*. 

For m3 = 1 we have 9* = sin^^(-\/2 — 1). Evaluating these constants numerically 
gives v{0) < —1.3 < as required. 

The branch 7 requires more work. First, it must be shown that 7 remains in the 
w > part of the energy manifold at least until it crosses the line 9 — 0. Thus is 
equivalent to v{9) being defined for 9 G [9* — tt,0]. Then we also require that that 
vi — w(— f ) < and = w(0) > 0. Recall that solutions can only leave {w > 0} in 
region Rjj. If v{9) is defined for 9 G [9* - tt, -|] and if ^^(-f ) < then in Rjj, 7 is 
trapped between two branches of stable manifolds and cannot reach the boundary 
curve {w ~ 0} (see figurejs]). Thus it suffices to show that f(— f ) < and v{0) > 0. 

To show that f(— f ) < 0. Recall that the initial value is v{9* — vr) = ~v* 
where w* = ^/6 « 2.44949. First we show v{-^) < -1.6 then we show that the 
change of v{9) on [— "^j- f] is at most 1.56. For the first part, use the estimate 
V{9) < V{-^) < 4 for 6* e [9* - tt,-^]. From ^ we have 



dv ^ 2|cos6i|V8- w2 
d9 - 1 + sin^ 9 



If w(— ^) > —1.6 we would have 



Sir 

dv f~ 4 2|cos6l|d6' 

< 



V8 - Je'-TT l + sin2 6' 
Both sides can be integrated exactly and the resulting formulas evaluated numeri- 

StT TT 1 

4 ' 2J 



cally to arrive at the contradiction 0.4459 < 0.4456. For 9 € [-^, - f ] we use the 



crude approximation 



dv 2^2W{9) 
d9 - 1 ' --"^ 



- sm 

Symbolic differentiation of the right-side shows that it is strictly decreasing on 
[— TT, — ^] which implies that its integral can be easily bounded using upper and 
lower Riemann sums. In particular, the upper sum for [— ^,— f] with 100 steps 
gives an upper bound of 1.559 < 1.56 for the change of v over this interval. 
Finally, we show that v{0) > 0. From (48 1 we have 

dv 2\ cos9\\/v*'^ — 
d9 - 1 + sin^ 9 ■ 

This can be used to show that «(— f) > If this were false we would have 

75 dv 2\cos9\d9 
^v*^ - i;2 ^ 1 + sin^ 9 ' 

The left side is ^ and using sin(6'*) = — sin(7r — 9*) = \/2 — 1 one finds that the 
right side is also ^, a contradiction. Therefore, t;(|) > — 
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Since v{6) is increasing we have v{9) > w(— f ) > = for 9 e [—-1,0]. 

If we had v{0) < then we would have v{t)'^ < 3 for 6* e [-|, 0]. Then we would 
have 

dv ^ 2\cose\y/v*^ -v{ty ^ 2|cos6l|\/3 

d0 ~ 1 + sin^ 9 ~ 1 + sin^ 6 ' 

Integration shows that v{d) increases by at least f\/3 for 6 E [— 1,0] giving the 
estimate 

w(0) > v{-l) + f \/3 > (f - l)VS > 
a contradiction to the assumption that v{0) < 0. Hence ?;(0) > as required. 
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